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Abstract: Hofava's proposal for non-relativistic quantum gravity introduces a preferred 
time foliation of space-time which violates the local Lorentz invariance. The foliation is 
encoded in a dynamical scalar field which we call 'khronon'. The dynamics of the khronon 
field is sensitive to the symmetries and other details of the particular implementations 
of the proposal. In this paper we examine several consistency issues present in three 
non-relativistic gravity theories: Hofava's projectable theory, the healthy non-projectable 
extension, and a new extension related to ghost condensation. We find that the only 
model which is free from instabilities and strong coupling is the non-projectable one. We 
elaborate on the phenomenology of the latter model including a discussion of the couplings 
of the khronon to matter. In particular, we obtain the parameters of the post-Newtonian 
expansion in this model and show that they are compatible with current observations. 
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1. Introduction 

Recently P. Hofava argued that it may be possible to construct a consistent renormalizable 
theory of quantum gravity within the framework of perturbative quantum field theory 
(QFT) Q (see also [g]). The proposal exploits the improved ultraviolet (UV) behavior of 
non-relativistic QFTs possessing an UV fixed point with anisotropic scaling of space and 
time. In this type of theories the UV behavior of the field propagators is improved thanks 
to terms with higher spatial derivatives. At the same time the number of time derivatives 



in the Lagrangian remains equal to two allowing to bypass the problems with ghosts arising 
in Lorentz invariant higher-derivative theories of gravity |J. Evidently, by breaking the 
symmetry between space and time one sacrifices Lorentz invariance (LI). The latter is no 
longer a fundamental property of the theory and may only emerge at low energies as an 
approximate symmetry. For matter fields in flat space this does not pose an immediate 
problem, other than the stringent observational constraint that Lorentz invariance needs 
to emerge to an extremely high accuracy, see |I], ||, ^, ^, [8|, |9| and references therein for the 
discussion of experimental bounds on Lorentz violating extensions of the Standard Model. 

However, abandoning relativistic invariance has dramatic effects for gravity. The rea- 
son is that in general relativity (GR) Lorentz symmetry is part of the local gauge group 
whose role is to remove the unphysical degrees of freedom contained in the metric, leaving 
only two massless helicity-2 modes for the graviton. Thus one expects that abandoning LI 
will lead to the appearance of new degrees of freedom in the theory. As already pointed 
out in Q this is indeed the case: in addition to the helicity-2 modes the Hofava model 
propagates a helicity-0 excitation 1 . The physical meaning of the new mode can be under- 
stood as follows. From the geometrical point of view the introduction of the preferred time 
coordinate amounts to equipping the space-time manifold with a foliation by space-like 
surfaces. In the gravitational theory the foliation inevitably becomes dynamical together 
with the geometry of the manifold. The helicity-0 mode is nothing but the excitation of 
this foliation structure. In this sense, the extra mode describes fluctuations of the global 
time, so we coin for it the name 'khronon' 2 . 

Importantly, the new mode does not have a mass gap and thus cannot be consistently 



decoupled at low energies [11, 12 1. This implies that, contrary to the original expectation 
HJ, the theory cannot flow to GR at low energies. Instead, one may entertain the possibility 
that at low energies the theory reduces to a (Lorentz-violating) model of modified gravity, 
with the modifications being small enough not to contradict the experimental data. The 
studies of modified gravity models, both in Lorentz-invariant jl4|, 15, [l(], 17] and Lorentz- 



violating [18, 19, [2(]] contexts, have shown that the properties of the extra degrees of 
freedom can make them fail as phenomenologically acceptable effective field theories (EFT) . 
This is precisely what happens in the original realization m of Hofava's proposal where 



the behavior of the khronon turns out to be pathological [21, g^, 12, £3], 24]. 

It is worth clarifying whether these pathologies are completely generic and invalidate 
the approach to quantum gravity proposed in ]!]], or they represent merely a failure of 
a specific realization of the general framework. The purpose of the present paper is to 
address this question. 

To this end we consider three different models of non-relativistic gravity. First we 
will reexamine the 'projectable' version of the original proposal [|| outlining its problems 
associated with the scalar sector. Then we will consider the two extensions of the original 



1 The helicity-0 mode is absent in a special case when the model obeys anisotropic Weyl invariance. 
However, this symmetry is not compatible with observations. For another proposal to eliminate the helicity- 
mode see lid]. 

2 From Greek xP 0V0( i ~ time. This should not be confused with the term 'chronon' appearing in the 
sense of "fundamental interval of time" in some theories |13l . 



proposal suggested in 12] to remedy these problems. The first one is based on a smaller 



symmetry group and can be viewed as a (power-counting) renormalizable version of ghost 



condensation [18]. Such a model allows to alleviate the problems of the low-energy EFT. 
However, we find that this leads to other troubles. The reason is the unavoidable presence 
in this formulation of a second helicity-0 mode at high energies that leads either to fast 
instabilities or to the break-down of the perturbative description. Hence, we conclude that 
this attempt does not offer a promising candidate theory of quantum gravity. 

The second extension is the one presented in pq| . In contrast with the two previous 
models, in this case the scalar sector is free from pathologies. A preliminary study of 
this healthy extension was reported in [^, |26|] . Here we analyze it in more detail and 
confirm that the model is compatible with phenomenological constraints for suitable choices 
of parameters. This suggests that it may serve as a starting point for constructing a 
renormalizable theory of quantum gravity. 

All the models we will consider contain flat space-time as a consistent background and 
are naively power-counting renormalizable. However, the behavior of the extra mode(s) 
is drastically different in the three cases, the differences stemming from the symmetries 
or other details of the particular realization. Hence, before going to our analysis let us 
spend a few words concerning the main features of the different possible implementations 
of non-relativistic quantum gravity. One can distinguish the following: 

• The choice of the anisotropic scaling to implement power-counting renormalizability. 
Namely, one postulates the scaling transformations 

xi-)-6~ 1 x, t^tb'H, (1.1) 

with a given critical exponent z, together with the scaling weights of the different 
fields. Then one classifies the operators in the theory according to their dimensions 
with respect to this scaling. The theory is (power-counting) renormalizable if the ac- 
tion contains only a finite number of terms of zero (marginal operators) and negative 
(relevant operators) scaling dimensions 3 . The case of a relativistic QFT corresponds 
to z = 1. In this paper we stick to the choice z = 3. As discussed in |Q|, this is 
the minimal value of z that allows to construct a power-counting renormalizable La- 
grangian for gravity in (3 + 1) dimensions. Larger values of z are also possible and 
lead to super-renormalizable models. 

• The subgroup of the four-dimensional diffeomorphisms (Diff ) under which the theory 
is invariant. The distinction between space and time enforces a preferred frame, and 
so a preferred time coordinate. As stated above, this corresponds to endowing the 
space-time manifold with an additional structure: a preferred foliation by space-like 
surfaces. In particular this means that the arbitrary reparameterizations of time 
t i—t- t(t,x) are not an invariance of the theory. Instead, the following unbroken 
symmetries have been considered in the literature: 



3 The dimension of the spatial coordinates xt is taken to be — 1. 



(i) xH>x(t,x) and t*-+i(t). (1.2) 

We will refer to these transformations as 'foliation-preserving Diffs', or FDiffs for 
short. This is the largest possible unbroken gauge group. It is the one originally 
considered in |ij. 

(ii) x h> x((, x) and t *-+ t = t + const. (1-3) 

We shall refer to this as the 'restricted-foliation-preserving Diffs' (RFDiffs). This 
symmetry arises in a number of effective field theories of modified gravity, such 
as the shift-symmetric A:— essence [27] or the ghost condensation [pLSjj . The in- 
variance under time translations implies existence of a conserved energy. 

(iii) x i—)- x(£,x). 

This is the unbroken group in potential-driven inflation and in non-shift-sym- 
metric £;— inflation around time-dependent spatially homogeneous solutions. It 



serves as the basis of the effective field theory of inflation [28, 29]. The action in 
this case contains explicit time dependence and there is no energy conservation. 

(iv) x I—)- x(x) and 1 1— >■ t = t + const. 

This option is realized in Einstein-aether theory pOf (see pl[ for recent review) 
and gauged ghost condensation [Q. It leads to the appearance of propagating 
helicity-1 degrees of freedom (in addition to the helicity-2 and helicity-0 modes). 

(v) x \— > x = x + £(£) and 1 1— > t = t + const. 

Here £(£) is and arbitrary time-dependent three-vector. This is the symmetry 



group of Lorentz violating massive gravity [19, 33 1. 



In principle it may be possible to construct power-counting renormalizable theories of 
gravity with any of the above unbroken symmetries. A natural expectation is that the 
larger the unbroken gauge group, the more constrained the model is, and the fewer 
degrees of freedom it contains. Investigating all the possibilities is beyond the scope 
of this article. We shall limit the analysis to the FDiff-invariant and RFDiff-invariant 
theories. In these cases there are only extra scalar modes, but their number grows 
when relaxing the symmetry. We will see that the covariant form of these theories 
involves only a single scalar field (a khronon) in addition to the metric. 4 It seems 
therefore appropriate to refer to these as "khrono-metric" theories. 

Finally, one may impose additional restrictions on the action that do not follow from 
a symmetry. Examples of such restrictions are the 'projectability' and the 'detailed 
balance' conditions of the original paper [J1J. In the present work we do not impose 
the detailed balance condition, and we consider models both with and without the 
projectability property. 



4 The equation of motion for the khronon field turns out to be of different order in time-derivatives 
depending on the symmetry. Hence the different number of propagating modes in the FDiff- and RFDiff- 
invariant cases. 



The paper is organized as follows. In Sec. |2|we introduce the basic notations and tools 
for our analysis. In Sec. |3]we consider the 'projectable' version of the original Hofava model 
and discuss problems associated with the extra scalar mode. In Sec. ||| we elaborate on an 
attempt to fix these problems by relaxing the symmetry from FDiffs down to RFDiffs. We 
show that the model obtained in this way does not provide a consistent candidate for quan- 
tum gravity. The relation between this model and the ghost condensation ][0| is discussed. 
In Sec. H we return to the case of full FDiff-invariance, now without the projectability 
condition. In this case the original action of Hofava's model must be supplemented by 



additional terms 25] allowed by the symmetry and power-counting renormalizability. In 
Sec. [0| we demonstrate that these terms make the scalar graviton stable and weakly cou- 
pled. In the rest of Sec. |5| we analyze the phenomenological bounds on the healthy model. 
In doing this we exploit the analogy [^, [34|] between the low-energy limit of the model and 
the Einstein-aether theory [[3(], 31]. In Sec. |6]we summarize our results and discuss future 



directions. Some details of the analysis are deferred to the Appendices. 

Readers interested in the phenomenology of the healthy model may skip Sees, y, |] and 
go directly from Sec. § to Sec. ||. 

2. General Setup 

2.1 Three theories under scrutiny 

In this section we introduce the basic ingredients of the Hof ava-type theories that we are 
going to consider. The field content includes the spatial metric jij, the shift N{ and the 
lapse N entering into the (3+1) (ADM) decomposition of the 4-dimensional metric, 

ds 2 = {N 2 - NiN^dt 2 - 2N i dx i dt - - fij dx i dx j . 

These fields transform in the standard way under the 4-dimensional coordinate transfor- 
mations. 

Model I: Hofava's projectable FDiff gravity. 

We first consider the FDiff-invariant case. Let us also impose the additional requirement 
that the lapse is 'projectable', i.e. that it does not depend on space coordinates, 

N = N(t). (2.1) 

Note that this restriction is compatible with the transformation rules for the lapse under 
FDiffs, 

N i-> N = N-L . 
dt 

One writes down the following action with two time derivatives jj]], 

m 2 r 

Si = -^- J d 3 xdt^N (KijK lJ - XK 2 - Vj) . (2.2) 



Here Mp is the Planck mass; A is a dimensionless constant; Kij is the extrinsic curvature 
tensor for the surfaces of constant time, 5 

Kij = ^ (iij ~ ViNj - VjiVO , (2.3) 

K is its trace, K = Kijrf 3 ; and the 'potential' Vj depends on the spatial metric 7^ via the 
3-dimensional Ricci tensor Rij and its spatial covariant derivatives. One notices that the 
first two terms in ( |2.2D , which comprise the kinetic part of the action, are invariant under 
the scaling (|l . l|) with z = 3 provided 7^, iVj and A?~ scale as 



lij ^ 7ij 



AT,- h^ 6 2 iV,- , AT ^ AT . 



In other words, the kinetic terms are marginal with respect to the anisotropic scaling with 
z = 3. A power-counting renormalizable theory is obtained by including in the action all 
possible marginal and relevant terms. This corresponds to considering the most general 
potential Vj containing local operators with the scaling dimensions up to 6. There is a 



finite number of these terms, whose classification was performed in 11]. For our purposes 
it suffices to write schematically, 6 



Vj = -£R + M~ 2 {A X R 2 + A 2 RijIP + . . .) 

(2.4) 



+ M~ 4 (B 1 RAR + B 2 R ij R jk Ri + . . .) 



where £, ^4 n , -B n are dimensionless coupling constants and dots stand for all the possible 
inequivalent operators of the given dimension. Note that we have introduced here the mass 
M* for the scale suppressing higher-order operators; this scale may or may not coincide 
with Mp. In what follows we set the parameter £ to 1 which can always be achieved by a 
constant rescaling of the time coordinate. In the absence of matter this does not affect the 



physical content of the model. We will return to the general case £ 7^ 1 in Sec. 5.4 . 

Omitting higher-derivative terms in the potential and setting A = 1 one formally 
recovers the action of GR O. 7 However, as we discuss below the limit A — > 1, A n ,B n — > 
is discontinuous: instead of flowing to GR the model ( |2.2[) , (|2.4|) becomes pathological in 
this limit. This model will be analyzed in Sec. y. 

Model II: non-pro jectable FDifT gravity 



Let us now relax the projectability condition ( |2.ip , so that the lapse N is allowed to depend 
both on time and space. In this case, the object 

a z = N-%N , 



5 Throughout the paper we use lower-case Latin letters i,j,... for 3-dimensional indices. They are 
raised and lowered using the spatial metrics jij, 7 1 -'. The covariant derivatives carrying these indices are 
understood accordingly. 

6 Recall that in three dimensions the Riemann tensor is completely determined in terms of the Ricci 
tensor. 



7 Still, the resulting theory is different from GR because of the projectability restriction (2.1) 



transforms covariantly under FDiffs and has scaling dimension 1. To obtain a (power- 
counting) renormalizable theory one should allow the potential to depend on etj. Again, 
this dependence must be limited to local operators of dimension up to 6. Note that a^ can 
be consistently excluded from the kinetic part of the Lagrangian. Indeed, recall that each 
time derivative raises the dimension of an operator by 3. Thus operators of dimensions less 
or equal 6 can contain together with a, at most one time derivative. It is straightforward to 
check that, up to integration by parts in the action, there are three possible combinations 
(all having dimension 5): 



K ij a iaj , K ij ViQj , KV, 



a 



These operators are odd under T and CPT transformations. While in general it might be 
interesting to study the effect of such operators, in the present paper we forbid them by 
assuming the CPT invariance. Thus we arrive to the following action 

M 2 r 
S n = ^ d 3 x dt^N (KijRV - \K 2 - V n ) , (2.5) 



where 



V H =Vj-a cue? + M-^CiOiAa* + C 2 {a t a 1 ) 2 + C 3 a iaj R lj + . . .) (2.6) 

+ AC 4 (I>iaiAV + D 2 ( ai a 1 ) 3 + D z a i a i a i a k B? 1i + ...). 

Note that due to the positive scaling dimension of Oj the potential again contains only finite 
number of terms. The projectable model can be recovered from ( p.6| ) by taking the limit 
a — > oo which enforces the spatial gradient of the lapse to vanish, diN = 0. The potential 
(|2.6| ) was first proposed in |P5[ . The model with the action Q2.5|), ( |2.6[) will be studied in 
Sec. |. 

Model III: RFDifF gravity 

The third model we are going to consider corresponds to relaxing the symmetry from FDiffs 
to RFDiffs. This has a dramatic effect on the possible structure of interactions in the theory. 
The reason is that the lapse N is now a scalar under the symmetry group 8 . At the same time 
its scaling dimension is 0. Therefore, all dimensionless couplings in the Lagrangian may 
acquire arbitrary dependence on it without spoiling the power-counting renormalizability. 
Besides, the reduced symmetry allows and power-counting renormalizability requires to 
include into the action a kinetic term for N. Thus the most general action reads, 

Sin = ^Jd 3 xdt^fN ( y X 1 (N)(N-N i d i N) 2 + X 2 (N)K i:J K^-X(N)K 2 -V III ) (2.7) 

s Note that the reduced gauge invariance is compatible with restricting the lapse TV to a fixed value, say 
N — 1. Then iV drops out of the action. The only difference of the resulting theory from the projectable 
model is the absence of the integral Hamiltonian constraint, so locally the two theories are equivalent (cf. 



the discussion in the beginning of Sec. 3.2). Thus all results about local dynamics of the projectable model 
apply to this case. 



with 

Vin =V(N) - £(N)R - a(N) a { a l (2.8) 

+ M- 2 (Ax(N)R 2 + A 2 (N)R ij R i ^ + C^mAa* + C 2 (N)(a i a i ) 2 + . . . ) 
+ MZ A (B X (N)RAR + B 2 {N)R ij W k R{ + D 1 (N)a i A 2 a i + D 2 (N)(a i a i ) 3 + ...). 

Note that the form of the kinetic term for N is fixed by the invariance under RFDiffs ( p..3| ) . 
Besides, relaxing the symmetry down to RFDiffs allows to include the standard potential 
terms V(N) for the lapse. The presence of arbitrary functions of N in the Lagrangian 
makes this theory perhaps less attractive, since it involves an infinite number of coupling 
constants. In spite of this, we find it instructive to explore what happens in this type of 
extension. We study this model in Sec. Q. 

Let us end this section by introducing a few notations that will be common in the 
analysis of all the above models. We are going to study the dynamics of small metric 
perturbations in these models above flat background. Since the most worrisome modes are 
the helicity-0 excitations we concentrate on scalar perturbations of the metric which we 
parameterize as follows: 

N = l + <f>, (2.9a) 

Ni = -?Lb , (2.9b) 



/A 
Hi = <% - 2 fa - t*&\ V - 2 & -^E , (2.9c) 

where A is the spatial Laplacian. Finally, the following mass scales will appear in the 
analysis: 

Ml = (A - 1)M P , Ml = aM 2 P , M\ x = \ X M 2 P . (2.10) 

2.2 Stiickelberg formalism and the khronon field 

In our study we will make extensive use of the Stiickelberg formalism. This formalism 
allows to single out explicitly the extra degrees of freedom appearing due to breaking of 
Diff invariance and study their dynamics in a transparent way. In the present context, it 
amounts to rewriting the action of the theory in a generally covariant form at the expense 
of introducing a compensator field which transforms non-homogeneously under the broken 
part of the 4-dimensional Diffs. It is worth stressing that this procedure does not introduce 
new degrees of freedom into the theory, but just makes explicit those already present. See 



1 15, 19, 33 1 for previous applications of the Stiickelberg formalism in modified gravity 



theories. 

In the present paper we follow the approach of [12 1 . One encodes the foliation structure 
in a scalar field (f(x). Namely, the foliation surfaces are identified with the surfaces of 
constant tp, 

ip = const . (2-11) 



Clearly, to define a regular foliation the field ip must possess non-zero time- like gradient. 
In this sense <p defines an absolute time, hence we call this field 'khronon'. This identifica- 
tion allows a covariant definition of the foliation structure. The actions ( p.2[ ), ( |2.5| ), ( [2.7] ) 
correspond to the frame where the coordinate time coincides with ip, 

<p = t. (2.12) 

We will refer to this choice of coordinates as the 'unitary gauge'. 

The action for p in a generic frame is obtained by realizing that the objects appearing in 
fl2.2j ), (^1]|), ( ^7| ) are the standard geometrical quantities that characterize the embedding 
of the hypersurfaces (|2.11| ) into space-time. The central role in the construction of these 
quantities is played by the unit normal vector 9 

- -^ (2.13) 



.'X 

where 

X = g^d^d u <p. (2.14) 

Note that u^ is invariant under the reparameterizations of ip, 

<p»<p = f(<p), (2.15) 

where / is an arbitrary monotonic function. This reflects the invariance of the foliation 
structure under reparameterizations of <p. The time-dependent VEV of the khronon field 
ip breaks the product of the reparameterizations ( 2.151 ) and general covariance down to 



the diagonal subgroup which in the unitary gauge coincides with the FDiffs (|Tj). One 
concludes that the covariant form of a FDiff-symmetric theory must be invariant under 



( 2.15 ). On the other hand, in the case of RFDiffs the symmetry of the khronon action is 



reduced to the shift symmetry 

ip i— > (p = ip + const , (2-16) 

allowing general dependence of the action on the derivatives of ip. 

Other geometrical quantities associated to the foliation are constructed out of u^ and 
its derivatives. We have the following expressions for the spatial projector: 

the extrinsic curvature: 

K.^ = P PP S 9 u v = -LpPPZV p V*<p , (2.17) 

VA 



the acceleration of the congruence of curves normal to the foliation: 

a M = vTVvUn (2.18) 



9 The Greek indices fi,i/,... are raised and lowered using the 4-dimensional metric g^ v . The same 
correspondence applies to the covariant derivatives carrying these indices. 



and the intrinsic Riemann tensor: 

<7?M _ pv pP pi pS (4) pa , irpir _ ir^ir 

,x - vpcr — r a r v r p r a ■"■ (3-y& ' '^p'^-ua l^-cr^vp > 

where in the last equation ^R a g„g is the 4-dimensional Riemann tensor. Now it is straight- 
forward to obtain the covariant form of the actions ([T^), Q2.5| ), ([2.7D by identifying the 
quantities appearing in the ADM decomposition with the appropriate combinations of d^ip, 
u p, Ppv, K-pv, etc. Namely, writing the covariant objects in the unitary gauge we obtain: 



X = N~ 2 , 




u = N , Ui = , 




pOO _ pOi _ q pij _ 


- 7 *J 


/V'2 7 — -**■ ij j C ^ J ^- / * 





In this way one ends up with an action describing GR plus a derivatively coupled khronon 
field (p. 

Postponing the detailed analysis to the following sections, let us discuss here an im- 
portant issue about the counting of degrees of freedom in the Stiickelberg picture. Note 
that the higher spatial derivatives of Hofava theories in the unitary gauge translate in 
the covariant picture into higher covariant derivatives of the khronon field. Consider, for 
example, the covariant form of the action ( |2.5D 

8„ = & 



f d 4 x^/^{ - WR - (A - 1)/C 2 - a a^ + ...}, (2.19) 



where /C = /C^, and K.^ v and a M depend on the khronon field via the expressions ( 2.17] ), 



(2.18|), (p. 13); dots stand for the higher-order terms 10 . It is easy to see that the second and 



the third terms in (2.19) contain four derivatives 11 of (p. Thus in a general coordinate frame 



the equations of motion for <p are fourth-order in time. Naively, this implies presence of 
ghosts and hence inconsistency of the theory. However, as discussed in |1^] this conclusion 
would be incorrect due to the following reason. Though the action Q2.1S ) has a covariant 



form, the physical content of the model does depend on the choice of the coordinate frame. 
The counting of degrees of freedom of the model must be performed in the preferred frame 
which is set by the preferred foliation. In this frame the number of time derivatives in the 
equations following from ( 2.19J ) is reduced to two, and the ghost modes are absent. 



While the above statement is obvious in the unitary gauge where perturbations of the 
khronon field are set to zero, it is rather non-trivial for general gauges where (p is allowed 
to fluctuate. Let us formulate this statement more precisely. Consider perturbations of the 
field p around a background (p, 

<P = V + X ■ 

We do not make any assumptions about the background (p, in particular, it does not need 
to obey equations of motion. For the sake of the argument we treat the metric as external 

10 We remind that we have set the constant £ in the potential (2.4) to 1. 
lr The higher order terms omitted in (2.19) contribute with even more derivatives. 



10 



assuming that it remains unperturbed; again, the background value g^ u (x) can be arbitrary. 
Let us now choose a coordinate frame such that the time coordinate coincides with the 
background value of the khronon field, 

t = ip . (2.20) 

Then in this frame the linearized equations of motion for \ following from (2.19) are second 



order in time. 

This proposition is proved in Appendix [A]. As explained there, the proof uses in an 
essential way the invariance of the action under the transformations ( p. 15 ). As a conse- 



quence it also applies to the projectable model Q2.2| ) which, being FDiff-invariant, obeys 
this symmetry. On the other hand, the above statement does not hold for the model ( |2.7| ) 
whose covariant form is not invariant under (|2.15). We will see in Sec. that this makes 



the equation of motion for the khronon in the model ( [2.7D fourth-order in time; as one can 
anticipate, this will lead to certain pathologies. 

3. Hofava's projectable model 
3.1 Gradient instability 



In this section we will study the projectable case (2.2), where the lapse N is assumed to 



obey the condition (2.1). Let us analyze the dynamics of small perturbations around flat 
background. At this stage it is convenient to work in the unitary gauge. We concentrate 
on the scalar perturbations of the metric and use the decomposition (|2.9| ), where due to 
the project ability condition <j> depends only on time. Substituting these expression into the 
action (|2.2j) we obtain the following quadratic Lagrangian: 



Cf ] = — £ j - 2ip 2 - 2V>A?/> + 4i/j\/AB + 4ipE 



(A - 1) (VAB + E + 2^) 2 - ^(AVO 2 - 1§4 V>A 3 V> 



Mr Mi 



(3.1) 



where the constants f\ and g\ are related to the coefficients in the potential (2.4); the precise 
form of this relation is not important for us. Note that the perturbation of the lapse (j) 
drops out of the quadratic Lagrangian because of the projectability condition: when <f> does 
not depend on the space coordinates all terms containing it are total derivatives. 
Integrating out the non-dynamical fields B and E, one finds 

Let us discuss this result. The positivity of the kinetic term for the field ip, needed for the 
positivity of kinetic energy, imposes the restriction: 

^i>0. (3.3, 



11 



On the other hand, the dispersion relation following from (|3.2| ) has the form 



2 _ ^ ~ 1 ( 2 h 4 , 91 6 



^ = -—M-m»^m*\- < 3 - 4 > 



We see that for A satisfying ( |3.3| ) the field ip is tachyonic at small values of spatial momenta. 
This instability can be cut off by the second term in the bracket at p ~ M*. One estimates 
the maximal rate of instability as 

Imw ~ a/I a -!| m * • (3-5) 

To obtain an acceptable phenomenology, we must require that the instability does not 
develop during the life-time of the Universe, 12 

V|A-1|M» < # , (3.6) 

where Hq is the present Hubble expansion rate. Recall that the physical meaning of M* is 
that of the scale suppressing higher-derivative operators in the gravitational action. In other 
words, it is the genuine scale of quantum gravity in the model at hand. A theoretically 
motivated value for M* would be a microscopic scale of order or somewhat below the 
Planck mass. However, to be as general as possible, we choose not to impose any a priori 
assumptions about the value of M*. Then we are left with the experimental lower bounds 
on M*. The direct bound is rather mild and comes from the tests of Newton's law at the 



distances ~ 10/xm [37| and reads, 

M* > 0.1 eV . (3.7) 



Even taking M* at the lower end of the allowed range, the stability condition ( |3.6| ) implies 
that the parameter A must be extremely close to 1, 

|A-1|<HT 61 . (3.8) 



This is unappealing from the theoretical viewpoint. In particular, it is unclear if 

can be preserved by radiative corrections. Most importantly, though, we are going to see 



that the condition (3.8) introduces a more important problem: the scalar mode becomes 



strongly coupled at unacceptably low energies. 

3.2 Strong coupling 

In this section we extend the study of the projectable model beyond the linear level and 
consider self-interaction of the scalar mode. As we will show, these considerations rule out 
the model as a weakly coupled alternative to GR. 

We will use the Stiickelberg formalism described in Sec. |2.2| . To construct the appro- 
priate khronon action one notices 12] that the projectable model (|2.2|) with the condition 



It was argued in [B5, 3q] that the instability may be cut off by non-linear effects which may lead to 



weakening of the bound (3.6). However, this scenario looks highly non-trivial and a detailed analysis of the 



non-linear dynamics of the model is needed to see if it is indeed realized. Thus we prefer to stick to the 



bound (3.6) following from the linear theory. 
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fl2.l| ) is locally equivalent to the theory with the same action where the lapse N is set 
to a constant; for example one can choose 13 N = 1. Indeed, the only difference between 
the two theories is in the equation obtained from the variation of the action with respect 
to the lapse. But this equation in the projectable model has the form of a Hamiltonian 
constraint integrated over the whole space and does not affect the local physics, cf. [38]. In 



Ref. [12] the condition N = 1 was implemented by introducing a Lagrange multiplier, and 
it was shown that the classical low-energy dynamics of the khronon field are equivalent to 
that of a pressureless fluid |jj|]. In the present paper we adopt a different approach more 
suitable for the analysis of the quantum properties of the system. Namely, we implement 
the constraint by adding to the action (|2.2|) a potential term 



f , M 4 / 1 Y 

S M = j d 3 xdt^N— h^-i) 



and considering the limit M — > oo. In this limit the potential becomes very steep and 
enforces N = 1. We will keep the coefficient M finite in the intermediate calculations and 
take the limit only at the last step. 

Now we are ready to write down the covariant form of the action for the projectable 



model. Using the prescription of Sec. 2.2 we obtain, 



4 71/T2 / X7U,^X7V,„ \ 2 



s, -\J**Fi \-m^r + ^x - 1). - f (r>- ^W) 

" (3.9) 
where ip is the khronon field, M\ and X are defined respectively in (|2.10|) and (2.14|). In 



this expression we have omitted the terms with higher spatial derivatives from the potential 
(|2.4| ). Note that fl3,9|) coincides with the action of the ghost condensate model |lq| . In 



particular, the last term in (3J3) is a covariant realization of the term with four spatial 
derivatives of the ghost condensate from Ref. [18] . However, whereas in [18] M and M\ 



are assumed to be of the same order, we are interested in the case M\ <C M. 

To get a first feeling about the khronon dynamics it is instructive to neglect the back- 
reaction of the khronon on the metric. Formally, this is achieved by considering the limit 

M P -> oo , M A , M - fixed . (3.10) 

The perturbations of the metric get frozen in this limit, so the khronon dynamics effectively 
decouples from that of gravity. Following the terminology adopted in massive gravity 
we refer to this regime as 'decoupling limit'. Note that the phenomenological constraint 
( |3.8j ) implies M\ <C Mp so that one of the decoupling limit assumptions is automatically 
satisfied. On the other hand the second assumption M <C Mp is in conflict with the fact 
that we want to take M eventually to infinity in order to recover the projectable model. 
This means that the proper analysis of the projectable model will require going beyond the 



13 This choice breaks the time reparameterization invariance. As already mentioned, an alternative way to 
implement the projectability condition, which preserves the full FDiff symmetry, is to consider the a — > oo 
limit of the non-projectable model. However, the approach adopted in this section allows to reach our goal 
in a more direct manner. 
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previous decoupling regime. However, working in the decoupling limit is still instructive. 
First, it will allow to outline the steps of the analysis in a simplified setting. Second, we 
will use the decoupling limit results in the next section where we extend the model to a 
version of the ghost condensation. With this in mind, let us proceed for the moment with 
the assumptions (3.10). 



Setting the metric perturbations to zero, hav = 0, and writing (p = t + x we obtain: 



Siy - I d 4 x 



M 



2 



M A X 2 - ^(A X ) 2 - M 4 x(diX? ~ Ml x((A X ) 2 + 2<9 a d i A X ) + . . . , 

(3.11) 
where we have written only the most relevant (cubic) interaction terms. This action de- 
scribes a non-relativistic scalar with derivative self-coupling. The dispersion relation for x 
reads, 

9 M l a , 

" 2 = ^ 4 • < 312 > 

We can compare this dispersion relation with that obtained previously for the scalar mode 
of the projectable model, Eq. (J3~f). In deriving ( t3.1l|) we neglected the higher-derivative 
terms. This explains the absence in ( |3.12| ) of terms with higher powers of momentum 



suppressed by M*. Moreover, in the projectable limit M — > oo (|3.12| ) reduces to u = 
which coincides with the | A — 1 1 — > limit of (|3.4|) . The latter corresponds to the decoupling 
limit flJTloD . 



From the presence of derivative interactions in ( [3.11 ) one expects the field x to become 



strongly coupled above certain energy/momentum scale. The easiest way to identify this 
scale is to perform power-counting. Due to the non-relativistic form of the dispersion 
relation for the x-held the power counting rules in our case are different from the standard 
ones. To identify the proper rules we follow the approach of [^, [T^] (see also |39(j). One 
notices that the quadratic part of the action is invariant under the scaling transformations 



x H- 6 _1 x , t H- b~ 2 t , x ^ b 1/2 



X 



The interaction terms in ( |3.11| ) scale as b 1 ' 2 and b 5 ' 2 implying that the strength of the 
interaction grows at short distances. To simplify the further analysis let us choose the 
units in such a way that the quadratic part of the action takes canonical form, i.e. it 
contains only order-one numerical coefficients. This is achieved by the rescaling 



M 2 M^H , X = M-Hl x 1/2 X 



which yields, 
S Ix = d 3 xdi 



X 2 ~ ( -^- ~ MM~ 3/2 mx) 2 - M- 3 M A 1/2 ^((Ax) 2 + 2d lX diAx) + . . 

Now the momentum and frequency cutoffs (strong coupling scales) of the theory are iden- 
tified as the appropriate powers of the scale suppressing the interaction terms, 

A p , dec = mm{M~ 2 Ml M^M~ 1/5 } , (3.13) 

A*, dec = min{M- 4 M A 6 , M 12 / 5 M" 2/5 } , 
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where the subscript "dec" reminds that these expressions are obtained in the decoupling 
limit. Here Cj is the frequency corresponding to the rescaled time t; going back to the 
physical frequency one obtains 



A,, 



dec 



mm 



{M- 6 M 7 X , M 2 / 5 M^ /5 } 



(3.14) 



Note that the cutoffs in momentum and energy are different reflecting the non-relativistic 
nature of the theory; they are related by the dispersion relation ( |3 .12 ). 

In the projectable limit M — > oo both cutoffs ( p.l3| ), (3.14) vanish which would imply 
that the model is strongly coupled at all scales. However, as already pointed above, this 
reasoning has a caveat: the expressions ( p.!3| ), (3.14) are obtained under the assumptions 
Mp — > oo (which in particular implies M <C Mp). We are going to see that a finite value of 
Mp raises the cutoff from zero; but it is still too low to be phenomenologically acceptable, 



cf. |23]. 

To obtain the correct value of the cutoff we have to take into account the fluctuations 
of the metric. Let us first see how this affects the analysis at the quadratic level. Taking 
the metric perturbations in the form ( |2.9| ) and expanding the action (|3.9| ) to quadratic 
order we obtain, 



S 



(2) 



d 4 x 



M 2 P 



( - 2ijj 2 - 2^A^ + 4c/>AV> + 4ipVAB + AipE) 



,2 M\ 



+M\x - 0) 2 - -A (VA£ + E + 2V + A X ) 



(3.15) 



Note that the first line here is nothing but the quadratic part of the standard Einstein- 
Hilbert action. If one fixes in ( 3.15 ) the gauge x = an d takes the limit M — > oo forcing <f> 
to be constant one recovers the unitary gauge Lagrangian 14 ( 3A ) . For our present purposes 
it is more convenient to choose instead the gauge 



B = , 2ip + E = 



where khronon perturbations do not vanish. To integrate out the non-propagating degrees 
of freedom recall that we are working in the regime M\ <C Mp. From the dispersion relation 



(3.4) we know that in this case the frequency of scalar perturbations is much smaller than 
the spatial momentum, 



"•• - Wp |P| « |P| • 



(3.16) 



Thus we can neglect the term tp 2 compared to ipAtp in ( 3.15 ). Then the equation of motion 
of tp implies tp = (p; substituting this into the action we obtain, 



S 



(2) 



<?X 



Ml 



+ M\ X ~ <t>? 



Ml 



(Ax) S 



14 This procedure does not reproduce the last two terms in (3.1). These higher-derivative terms have been 
omitted in the low-energy action ( 3.15). 
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Varying with respect to eft one finds 



M 4 x 



MlA + M 4 



and hence 



SJ 2) = / d'x 



MpM 4 
M 2 A + M 4 



X&X 



Ml 



(Ax) 5 



This gives the dispersion relation 



ur 



M A 2 
2A% P + 



2M 4 ' 



(3.17) 



(3.18) 



(3.19) 



which reduces to ( |3.12| ) in the decoupling limit Mp — > oo. On the other hand, in the 
projectable limit M — > oo it correctly reproduces the first term in the exact dispersion 
relation (|1|). 

Let us turn to the interactions. One has to compare contributions coming from various 
terms in the action (|3.9|). We start by considering the second term in (|3,9|). Due to the 
inequality ( 3.16| ) the leading interaction is given by the cubic term with the smallest number 
of time derivatives. A simple calculation yields 



S 



(3) 

M 



d 4 x[-M 4 ( X -Md iX ) 2 + ...] 



<fx 



MpM 4 

mJaTm 4 



&m x ? + ■■■ 



(3.20) 



where in the second line we substituted 4> from ( |3.17| ). Note that in the decoupling limit 
this interaction term reproduces the third term in the action ( |3.11| ). Next, let us estimate 
the contributions from the first and third terms in ( |3.9| ) and show that they are suppressed 
compared ( [3.20 ), This can be done without explicit calculation on purely dimensional 
grounds. Indeed, ( 'R does not explicitly depend on the khronon and has mass dimension 
two. Thus the possible leading contribution (the one without time derivatives when written 
in terms of the metric) has the schematic form 



Mp^ 



M]xj) 2 A(j) 



Mpx 2 A X 



where in the last relation we used Eq. ( |3.17j ) and took the limit M — > oo. This contains 
more time derivatives than ( J3.2CJ ) and hence is suppressed. Finally, the schematic form 
of the leading contribution from the third term in ( |3.9D is M|x(^x) 2 - This is clearly 
suppressed compared to ( 3.20|) by the ratio (M\/Mp) 2 . 

In the projectable limit M — > oo both the quadratic action (3.15) and the interaction 
term fl3.20|) simplify. Combining them together we obtain 



S, 



d 4 x 



M 2 P x&X 



Ml 



(A X ) 2 - Ml A x (d iX ) 2 + 



It is now straightforward to find the strong coupling scales of the model. Upon bringing 
the quadratic part of the action into canonical form by the rescaling 



t = MpM'H 



X 



Ml 



V2 Mr l/2 



X , 
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one finds that the interaction term is suppressed by the parameter 

A p = M A 3/2 M p 1/2 , 

which thus sets the momentum cutoff of the theory. The energy cutoff contains additional 
factor M\Mp and reads 

A,, = M 5 x /2 M p 3/2 . 

Note that these cutoffs go down to zero in the naive GR limit A — > 1 (M\ —> 0). This 
agrees with the results of [23|. 

Using the definition of M\, Eq. ( p.lOj) , we find that the stability bound ( |3.£| ) implies 

A p <l(T 17 eV. (3.21) 

This corresponds to the distance of order 10 13 cm, within which the theory is strongly 
coupled and the perturbative analysis breaks down. Note that the strong coupling cannot 
be resolved by the higher-derivative terms: such terms would lead to sizable modifications 
of the Newton's law that are forbidden at these distances. We conclude that the projectable 
model fails to provide adequate description of gravity within distances ~ 10 13 cm. 

4. Extension a la ghost condensation 

The analysis of the previous section suggests a possible way to address the problems of the 
projectable model. The idea is to consider a modification of the theory with finite value 



of the parameter M in the khronon action (|3.9|) [12]. This approach implies reducing the 
symmetry of the theory from FDiffs down to RFDiffs as the khronon action is no longer 
invariant under the general reparameterizations ( 2.15| ) but only under the shifts ( 2,16| ). 



We presently explore this option and show that in spite of the possibility to improve the 
behavior of the khronon at low energies the theory still encounters serious problems in the 
UV completion. 

As already pointed out, at finite values of M the low-energy khronon action fl3.S| ) 
coincides with the action of the ghost condensate model [18|. The latter is a consistent 



effective theory describing low-energy modification of GR. From the EFT point of view it 



is natural to assume that the scales M and M\ appearing in (3J)) are of the same order. 



Then both the momentum and frequency cutoffs obtained in the decoupling limit ( 3.13; ) 



( 3.14 ) are finite and of order M. From the dispersion relation ( j3.19|) one observes that the 



instability of the scalar graviton at low values of momenta is still present in this model but 
is truncated by the quartic term at 

p ~ M 2 /M P . 

This corresponds to the maximal instability rate 

Imw~ M 3 /Mp . 

This rate is slower than the present rate of Hubble expansion for 

M < 10 MeV, (4.1) 
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which would make the instability harmless. 15 

The problems appear, however, when one tries to UV complete the model to a renor- 
malizable theory of gravity. First, a necessary requirement to obtain a weakly coupled UV 
completion is that the strong coupling observed within low-energy EFT must be resolved 
by the higher-derivative operators, cf. [26|. This implies that the scale M* suppressing 



higher derivatives should be lower than M, M* < M. Thus we obtain 

M* < 10 MeV . (4.2) 

Note that this is many orders of magnitude smaller than Mp so that the model contains two 
very different scales. For the gravitational sector alone, this seems to be only an esthetic 



problem: the values fl4.2|) are compatible with the direct gravitational bound (3/7). Besides, 
we expect this hierarchy to be stable under radiative corrections. The scale which effectively 
truncates the power-law divergencies in the model is the lower scale M*; therefore, the 
corrections to both Mp and M* are small 16 . However, the real tension arises when we take 
into account coupling of the theory to matter. The size of the higher-derivative Lorentz 
violating terms in the action of the matter sector is experimentally constrained from below. 
The leading effect of these terms is the modification of dispersion relations of the matter 
fields at high energies due to contributions with higher powers of spatial momentum, see 
Eq. (5.2C) below. The observational lower bound on the scale M* a ' suppressing these 



contributions, Eq. ( 5.21 ), is much higher than 10 MeV. On general grounds one expects 



the scale M* a to enter into radiative corrections for M* thus giving rise to a fine-tuning 

problem 17 . A more detailed study of this issue is beyond the scope of the present article. 

The second problem stems from the reduction of the symmetry from FDiffs down 

to RFDiffs. To get a renormalizable theory one must consider the most general RFDiff- 

invariant action fl2.7| ) which contains, among other terms, the term with time derivatives 

of the lapse, 18 

M 2 r 

Sxl = ~t J d3x dt ^ N (Ar " NidiN ? ' ( 4 - 3 ) 

where M\ 1 is defined in (2.10). The easiest way to see that this term leads to pathologies 



is to adopt the Stiickelberg picture. The covariant form of (|4.3| ) reads 

4,4^25!!^:. ,4.4) 

Let us work in the decoupling limit; the conditions for the validity of this approximation 



will be specified below. One writes ip = t + x an d expands (4.4) up to quadratic order in x : 



Sai = ^ / d 4 x (x) 2 



15 It is argued in MOf that nonlinear dynamics of the ghost condensate suppresses the exponential growth 
of perturbations predicted by the linear theory, making the model phenomenologically viable for M up to 
100 GeV. 

16 See also the discussion of a similar issue in the next section. 

17 The fine tuning may be absent if there is a symmetry protecting M» from this kind of corrections. 

18 Even if this term is absent in the bare action, one expects it to be generated by quantum corrections. 
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This contains four time derivatives of the khronon perturbations which implies the presence 
of a second helicity-0 mode besides the excitation studied so far. As the new mode appears 
due to higher time derivatives one is tempted to qualify it as a ghost. In fact, we are going 
to show that this mode also exhibits a gradient instability, so it is more appropriate to call it 
tachyonic ghost. In Appendix [B] we discuss the difference between the tachyonic ghost and 
an ordinary ghost with stable dispersion relation 19 . In particular, we demonstrate there 
that the presence of the tachyonic ghost is compatible with the fact that in the unitary 
gauge the action for the model, Eq. fl2.7|) , is only second order in time derivatives. 

To proceed we need the quadratic action for the khronon perturbations. This is ob- 
tained in the usual way by first covariantizing the terms entering fl2.7| ), (|2.8| ) and then 
taking the decoupling limit. This yields, 



Silly 



(fx 



Ml 



4 -2 



Mxr + M 4 x 



Ml 



M 2 

(Axf + ^m) 



(4.5) 



In deriving this expression we have omitted the contributions of terms with higher spatial 



derivatives; thus the domain of validity of (4J3) is restricted to spatial momenta smaller 
than M*. The last term in (|4.5| ) represents the contribution of the term actio 1 in the 
potential fl2.8|) (recall the definition of M a in (|2.10|) ). The action ( [Q| ) describes two modes 
with dispersion relations 



u 



M 4 
Ml 



+ 



M 2 
2M 2 

Ai 



P 



NV 



/ M 4 



Mi 



+ 



M 2 
2M 2 

Ai 



P^ 



+ 



Ml , 

-irP 4 



M 



(4.6) 



Ai 



At small values of momenta we find, 



2 
u old 



oj7. 



Ml 



2M 4 

2M 4 



P^ + 



M 2 

Ai 



M 2 

a . 

M 2 ■ 

Ai 



2M 4P ^ 



(4.7) 



In the first expression one recognizes the dispersion relation ( ft. 12 ) for the previously en- 
countered gapless mode. The second mode (4/7) is unstable. It possesses a frequency gap 
of order M 2 /M\ 1 . For the choice M\ 1 ~ M, natural from the low-energy effective theory 
point of view, the frequency of this mode lies above the scale M. Thus this mode would be 
simply discarded as unphysical in the EFT considerations with cutoff M, such as in [18]. 



However, our case is different: we are looking for a UV-complete model and thus have to 
take into account all possible excitations of the system. 

According to the above expressions the instability rate of the new mode grows with 
momentum, the fastest instability occurring at |p| ~ M*. [At larger momenta terms with 
higher spatial derivatives in the full action ( |2.7| ) can, in principle, stabilize the mode.] A 
lower estimate for the instability rate is obtained by keeping only the last term under the 



square root in Eq. (4J5). This yields, 

Im uj. 



new 1-^j lvl * lvl \ 



1/2 



(4. 



9 Of course, both type of ghosts signal pathologies of the theory. 
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where we have assumed 

M~M a ~ M x > M* . (4.9) 

Requiring the instability to be slower than the Hubble rate Hq and taking into account the 



experimental bound (3.7) on M* one obtains, 

M Xl > 10 60 eV « W 32 M P . (4.10) 

Having such a huge value for M\ x is rather unsatisfactory. First, it means that one intro- 
duces into the theory, besides M* and Mp, one more hierarchical scale M\ x 3> Mp. We 
do not know whether this hierarchy is technically natural or not. More importantly, this 
term also gives rise to interactions which are enhanced by the large value of M\ t . We will 
demonstrate shortly that this reintroduces the strong coupling problem. 

Before addressing the interactions let us make a step back and discuss the conditions 
for the validity of the decoupling limit. The simplest requirement would be that all the 



parameters in the khronon action (4J3) are smaller than the Planck mass. This is automat- 



ically satisfied by the parameters M\, M, M a that must be small for the stability of the 



gapless (ghost-condensate) mode. However, the value ( 4.1C ) of M\ 1 clearly violates this 
condition. Thus additional considerations are needed to establish whether the decoupling 
limit holds or not for M\ 1 as large as ( |4 . 1 0| ) . A rigorous method to do this is to verify that 
the terms in the Lagrangian describing mixing between the khronon field and the metric 
components are small compared to the other terms. Instead of following this route we will 
take a shortcut and compare the decoupling limit dispersion relations (|4.6| ) for the scalar 
modes with the exact expressions; their coincidence will serve as a criterion for decoupling. 
The exact dispersion relations are obtained in Appendix [B] within the unitary gauge, see 
Eq. (|H|). In the case M Al > M P they simplify to 20 , 

J 1 - 

4Mp ^ V lQM p ' Af£ 




On the other hand, the decoupling limit result ([4.6D at large M Xl becomes 



.1 _ j. M A 2 



" = ±j£/ ■ <«*> 



The two expressions coincide and thus the decoupling holds, provided 



M Xl « ^ . (4.13) 



20 Note that these expressions immediately imply a lower estimate for the instability rate of the new mode 
which is complementary to (4.8), 



Imuj new > — — p . 

V2M P 

This estimate does not depend on M\ x . Evaluating it at jpj ~ M* and requiring the instability to be slower 

than the Universe expansion rate we obtain the upper bound 

A/ A <0.1eV. (4.11) 



This is marginally compatible with the requirement (4.9) and the experimental bound (3.7) 
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Note that this allows for M\ 1 to be much larger than Mp. Taking M\ ~ 0.1 eV as in 
the previous estimates (see also the footnote 20) one finds that the decoupling limit is 
applicable up to M\ 1 of order ( f4.1C ) . 

The previous considerations show that to estimate the scale of strong coupling we can 
concentrate on the self- interaction of the khronon neglecting its mixing with the metric. 
Clearly, for large M Al the leading interactions come from the term ( |4.4j ). At the cubic level 
we obtain the following contribution, 



S 



(3) 

Ai 



-2M( 



d x xdixdiX ■ 



This must be combined with the quadratic action ( fl.5| ) . The form of the dispersion relation 
(4.12) shows that at large M\ x the latter action is dominated by the first and the third 
terms. Retaining only these two terms we obtain the khronon action in the limit of interest 



Silly 



I 



d*x 



Ml 



(X? 






2Ml i xd iX d i x + ... 



The quadratic part is brought to canonical form by the rescaling 

t = M-^M^i, x = M^M-^ X . 

From the scale suppressing the interaction term in the resulting action one reads off the 
momentum and frequency cutoffs, 



A,. 



M A 5 / 4 M- 1/4 



A, 



m a 7/4 m a -; 3/4 



Substituting here M\ from ( [4.1 1| ) and M Al from ( |4.10[ ) one obtains 21 A p < 10 -16 eV. Such 
a low cutoff is phenomenologically unacceptable. This shows that the initial hope, namely 
that the model could be UV complete and weakly coupled, is not met in reality 22 . 

The general conclusion of this section is that reducing the symmetry of the theory from 
FDiff to RFDiffs cannot possibly improve the scalar sector. The reason is that the smaller 
symmetry allows for the new operators such as (^^) which bring in additional difficulties 
in the form of tachyonic ghosts. In retrospect, we learn that the symmetry under the 



FDiffs (or the khronon reparameterizations (2.15) in covariant language) plays the quite 
important role of preventing this kind of pathologies. In the next section, we return to the 
other model compatible with the FDiffs — the general non-projectable case — where this 
problem is automatically turned away. 23 



21 Note that this is close to the result in the projectable case, Eq. (3.21). 

22 One may wonder if it is possible to cure the model by choosing some of the parameters in the action 



(4.5) negative. A simple reasoning demonstrates that this is not the case. The analysis in the unitary gauge 
(Appendix 0) shows that negative values of MjJ , M\ imply negative energies in the UV. Next, making 
M 2 negative does not essentially improve the behavior of the tachyonic ghost. Finally, the choice AI 4 < 0, 
while removing the instability from the mode (|4.7|), destabilizes the old gapless mode. This brings us back 
to the situation studied in the previous section. 

23 Clearly, the symmetry (2.15) must be non-anomalous in order that this statement remains valid in the 



quantum theory. This question will not be addressed here. 
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5. The healthy extension 

5.1 Stability and absence of strong coupling 

We presently consider the non-projectable FDiff- invariant model with the action (|2.5| ), 
(|2.6|). We are going to show that, remarkably, the problems associated to the additional 



modes which plague the two models considered previously are absent in this case [25, p6[. 



Let us start by studying the spectrum of linear perturbations around flat background. 



As usual, we concentrate on the scalar sector and use the decomposition (2.9). In the 
unitary gauge the quadratic Lagrangian reads 



M2) Ml 

~ 2 



2V> 2 - 2?/>AV> + 4(/>A^ + AipVAB + AipE - (A - 1) (y/AB + E + 2^j 



2 



ft f\-i.\2 2 /2 A , A , h ,\i\2 91 i , A 3 I 2 52 , . ; 



(5-1) 

where f n , g n are related to the coefficients in front of the higher derivative operators in 
the potential fl2.6| ); the precise form of this relation is not important for our purposes. 
Integrating out the non-dynamical fields B, E and </> we obtain, 



, (2) __ M|f 2(3A-l) - 2 P[M~*A] 
2 I A-l w V Q[M- 2 A}' 



47 = -£\ \ i V + ^!,;- 2A i A^ \ , (5.2) 



where the polynomials P, Q have the form, 

P[x] =(92 ~ 9i9s)x 4 - [gih + 9sfi - 2£ 2 / 2 )x 3 + (/| - 4 52 - /1/3 - 2g 3 - gi a)x 2 

- (2/ 3 + fia + 4/ 2 )x + (4 - 2a) , (5.3) 

Q[x] =g 3 x 2 + f 3 x + a . (5.4) 



The Lagrangian fl5.2| ) describes a single mode which is free of pathologies provided that 
two conditions are satisfied. First, the positivity of the kinetic term can be achieved by 
choosing 24 A > 1. Second, from the dispersion relation of the propagating mode ?/>, 

W " 2(3A=I) Q[-P 2 /M2] P ' (5 - 5) 

one reads off the condition to avoid exponential instabilities, 

P[x]/Q[x] > at x < . (5.6) 

This puts certain restrictions on the coefficients a, f n , g n , that are presented in Appendix |C]. 
In particular, at low energies the dispersion relation takes the form, 



24 We do not consider the other option A < 1/3 because it corresponds to a strong deviation from GR, 
unacceptable from the phenomenological viewpoint. 
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Thus stability requires 25 

< a < 2 . (5.8) 



Note that the dispersion relation ( |5.7| ) describes a gapless mode propagating with constant 
velocity which is generically different from one (the velocity of the helicity-2 modes, i.e. 
gravitons). This signals that in the model at hand Lorentz symmetry is broken down 
to arbitrary low energies. In the UV the dispersion relation ( |5,5D takes the form lo 2 oc p 6 
which obeys the anisotropic scaling with z = 3. This is compatible with the power-counting 
arguments in favor of renormalizability. 

It is important to realize that the healthy behavior of the scalar mode can be achieved 
simultaneously with the stability in the sector of the helicity-2 perturbations. Consider 
operators in the action ( |2.5| ), ( ^^ ) which contribute at the quadratic level. Upon integrating 
by parts and using Bianchi identities one obtains a list of 10 inequivalent combinations, 

(dim 2) R, a-ia 1 , 

(dim 4) RijR ij , R 2 , RV t a\ a^Aa* , 

(dim 6) (V t R jk ) 2 , (ViR) 2 , ARV ia \ cjjAV . 

The dispersion relation of the helicity-2 modes depends only on the coefficients in front 
of the operators in the first column of this list. After fixing these to ensure stability of 
the helicity-2 modes, we still have the freedom to choose the coefficients of the remaining 
operators in the list. The number of free parameters matches with the number of coefficients 
a i fn, 9n in the scalar Lagrangian implying that we have freedom to adjust the latter 
coefficients to satisfy (|5.6|). 

To get more insight into the dynamics of the model (in particular, at the nonlinear 
level), we make use of the Stiickelberg formalism. The covariant form of the model action 
was given before, see (2.19). It is convenient to rewrite it as 



S H = —^ / d 4 xV=^{ (4) i? + (A - 1)(V M ^) 2 + a u»u u V^V u u p + ...}, (5.9) 



where u^ is related to the khronon field <p by (2.13|) and dots stand for the terms with 



higher derivatives which are not important at low energies. The action (|5.9| ) is closely 
related J2^, 34] to the action of the Einstein-aether theory |Kj (see Q] for a review). The 



latter is an effective theory describing breaking of Lorentz-invariance by a time-like vector 
field with unit norm. The difference between our case and the Einstein-aether theory is 
that in our model the unit vector u p is by construction hypersurface orthogonal, i.e. it 
is completely characterized by the scalar khronon field (p. The similarity between (|5.E| ) 
and the Einstein-aether theory will be exploited in the next subsections where we study 
phenomenological consequences of the model. 



Let us use (5.9) to study interactions of the khronon perturbations at low energies. 



It is convenient to introduce the scales M\ and M a as defined in ( |2.10| ). These scales 



25 As discussed in Sec. 2.1 



by taking the limit a — > oo in the model considered here one obtains the 
projectable Hofava's model. Clearly, the corresponding values of a are outside the stability range (p.q), so 
we again find that the projectable model is unstable. 
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characterize the khronon action. We will assume them to be much smaller than Mp, so 
that the metric perturbations are frozen out (in other words, we will be working in the 
decoupling limit). This assumption is justified by the phenomenological bounds ( 5.35| ) that 
will be obtained in Sec. 5.4, and which constrain the dimensionless parameters |A — 1 
to be much smaller than one 26 . Writing down <p 
order in x we obtain, 



o 



t + x an d expanding (5.9) up to cubic 



Sin 



d 4 x 



M 2 M 2 



+ 



M 2 xX ((&x) 2 + 2dadiAx 
M 2 (xdiXdiX ~ dixdjxdidjx 



+ 



(5.10) 



Let us analyze this expression. The action ( 5.10| ) describes a propagating mode with 
dispersion relation 

• 2 - x ~ 2 (5.11) 



0J 



Ml, 



This coincides with the exact dispersion relation (5/7) in the limit |A — 1|, a < 1. The 



form of the action (5.10) is uniquely fixed by the reparameterization symmetry (2.15) of 
the khronon and by the Lorentz symmetry, both non-linearly realized on the khronon 
perturbation x- Up to quadratic order the reparameterization transformations read 

X^X + £(t) + e(t) X + ... ■ 

where e{t) is an arbitrary function of time. Under boosts the field x transforms as 

X !->■ x + £i%i + ZiXiX + titdiX + • • • , 

where £« is a 3-dimensional vector characterizing the boost. It is straightforward to check 
that ( 5.10| ) is the only action invariant under these symmetries up to cubic order. Finally, 
one observes that ( [5.10D contains an interaction term with three time derivatives; it pro- 
duces a contribution with third time derivative in the equation of motion for x- One may 
be worried that this leads to appearance of a new unwanted degree of freedom. However, 
from the analysis in the unitary gauge, we know that this degree of freedom is spurious 
and it should be possible to eliminate it with an appropriate choice of variables (which 
corresponds to fixing the foliation consistently). To see how this is done explicitly, one 
considers the change of variable 

X = X + XX • 
This substitution can be interpreted as the change of the time foliation, t *-+ t — %■ Indeed, 
up to cubic order x(^ x ) = x(t ~ X> x )- I n terms of the new variable the action takes the 
form 



Siiy 



d"x 



^ m? ^A m? 



M 2 x Ax A X 



+ Mi 



^x(diX) 2 



dixdjxdidjx + 



(5.12) 



26 Note though that the bounds (5.35) are much weaker than the analogous bounds for the models of the 
previous sections, cf. (J3.8[) , (ft.l|). 
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We observe that a term with three time derivatives is still present but its structure has 
changed: its contribution to the equation of motion contains now only second time deriva- 
tive. One concludes that no new degrees of freedom appear. It is worth emphasizing that 
the existence of the change of variables with the above properties follows from the general 
statement proved in Appendix |A| that the equations of motion for the khronon field are 
second order in time for an appropriate choice of the coordinate system. 

The low-energy action ( |5.12 ) contains derivative couplings which become stronger as 



the energy / momentum grows. Let us estimate the scale where the strength of the inter- 
actions would become of order one. Performing the rescaling 

t = M a Ml x t, x = M-^M' 1 ' 2 x , 



which casts the quadratic part of the action into canonical form we obtain, 
S IIy = / d 3 x d£ 



mf my- M^. A . Ai 



M a 



3/2 



V2 ^1/2 

(5.13) 



M A 1/Z . . 2 M' a 



+ — \jKX(diX) TjKdixdjxdidjX + ■ ■ 

2Mi' M^ 

From the scales suppressing the interaction terms one reads out the momentum and fre- 
quency cutoffs of the low-energy description: 

A p = min {M~ l / 2 Ml /2 , M^M^ 12 } , (5.14) 

A w = min {M X J 2 M]I 2 , M" 3 / 2 M A 5/2 } . (5.15) 

If (|5.9| ) were the full action of the theory, it would become inconsistent at energies / 
momenta above these scales [JO]]. Note that A w and A p are related by the dispersion 
relation fl5.11|) and are different if the khronon velocity differs from 1. We see that the 
scales ( [5.14D , (|5.15| ) are low whenever there is a large discrepancy between M a and M\. 
On the other hand, in the case M a ~ M\ the cutoffs essentially coincide with M a , 

A p ~ A w ~ M a . 

We concentrate on this latter case in what follows. 

Of course, the action (|5.9[) , and hence (|5.13 ), represents only the low-energy part of 



the full action (2.5). So the existence of a finite cutoff for the low-energy theory (| 
does not imply any inconsistency. At energies/momenta larger than M* one has to take 
into account the higher-derivative terms in the full action, and any conclusions drawn from 



( |5.9| ) become invalid. By construction, the role of the higher-derivative terms in (2.5) is 
to modify the power-counting rules at high energies in such a way that all the interactions 
become marginal under the anisotropic scaling. One concludes that strong coupling is 
avoided provided the scale M* of higher-derivative operators is lower than 27 A PtU} [26]. 



In particular, one can show Pq] that in theories obeying the anisotropic scaling (1.1) with z — 3 in the 
UV the tree level unitarity, whose breaking usually serves as the signal of strong coupling, is automatically 
preserved provided A/* < A PjW . This is essentially due to the peculiar kinematics of these theories that 
makes the unitarity bounds milder at high energies as compared to the relativistic case. 
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This gives the upper bound 

M*<M a . (5.16) 

This is a remarkable inequality: it relates the scale of quantum gravity M* with the pa- 
rameter M a of the low-energy Lagrangian which can be probed experimentally. The exper- 
imental bound (|5.35| ) derived in Sec. |5.4| constrains M a to be lower than 10 15 -r 10 16 GeV. 
According to ( |5.16D this translates into the bound 

M* < 10 15 -T- 10 16 GeV . (5.17) 

Let us emphasize that, contrary to the claim in [41|, the choice of M* parametrically below 



Mp does not introduce a fine-tuning in the model. In fact, having M* somewhat below 
Mp is technically natural. From the point of view of the low-energy theory, the reason is 
that the cutoff for the power-law growing of the couplings is set by the lower scale M*. 
Thus neither Mp nor M* receive large quantum corrections. 

5.2 Phenomenological considerations 

The absence of notorious pathologies in the model (|2.5|) makes it worth to have a closer 



look on its phenomenological consequences. We are interested in the phenomenology of the 
model at energies much lower than M*. The natural language for this analysis is provided 
by the Stiickelberg formalism where at low energies the theory reduces to GR plus the 
khronon field. Then all non-trivial effects of the model are clearly interpreted as due to the 
presence of the khronon. Before proceeding we have to specify the coupling of the khronon 
to the fields of the Standard Model, to which we collectively refer as "matter" . 

The FDiff symmetry requires that matter couples to the khronon via geometrical ob- 
jects, such as Ufj,, a M , fC^, etc. Possible interactions fall into two classes having qualitatively 
different phenomenological consequences. The first class consists of couplings which contain 
derivatives of the vector u^; examples of this type of couplings are: 

a^ip , /C"V 7/ AV , (5.18) 

where ip is some fermionic matter field. Importantly, the combinations of the khronon 
field entering into the operators ( |5.18| ) have vanishing VEV in the flat background. There- 
fore, these operators start linear in the khronon perturbation \ an d describe its derivative 
interaction with matter. It is easy to see that in terms of canonically normalized fields 
these couplings have dimensions larger than four 28 and are suppressed by the high-energy 
scale M*. Given that M* is large, one expects the effect of these operators to be highly 
suppressed at the energies accessible to the present-day experiments. It is still possible 
that couplings of the type (|5.18j ) may be probed in some precision measurements. This 
issue is, however, beyond the scope of the present article. 

The second, more 'dangerous', class of operators describe coupling of the matter fields 
directly to the vector u M ; examples are: 

u^T'V » «*V<h/A</> , u^u^d^i) . (5.19) 



3 We use here the standard power-counting rules relevant for the low-energy physics. 
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The crucial property of these operators is that they give rise to Lorentz-violating effects 
within the Standard Model as they couple matter fields to the VEV of u«. Such effects 
and constraints on them have been extensively studied in the literature, see [||, ^, ||, || for 
reviews. The effects of the operators ( |5.19| ) fall into two categories: those that become 
pronounced only at high energies, and those that persist at all energy scales. 

Effects of the first category are associated to Lorentz-violating operators of dimensions 
larger than four such as the last operator in ( |5.19| ) . The strength of these effects is charac- 
terized by the scale suppressing the higher-order operators which we will collectively denote 
as M* a (one should keep in mind though that this scale can, in general, be different for 
different operators). In the model at hand it is natural to assume that M* a is of the 
same order of magnitude as the scale M* appearing in the gravitational sector. However, 
we stress that this is an additional assumption: in general the scales M* and M* a may 
be different, so we prefer to keep different notations for them. 

The experimental data constrain the scale M* a ' from below. A rather robust bound 
comes from astrophysical observations and exploits the fact that the higher-order operators 
lead to modification of the dispersion relations of the matter fields, in particular, photons, 
at high momenta 29 

D 4 

E 2 = m 2 + p 2 + r] 7 ,. + . . . , (5.20) 

> M {mat) 

where r\ is a dimensionless coefficient. This would produce a frequency dependent delay in 
the arrival times of 7-rays emitted by a distant source. The absence of such a time-lag in 



the signals coming from active galactic nuclei [44] and 7-ray bursts [45] yields the constraint 



(assuming the coefficient rj in ( |5.2C ) is of order one), 



M (mat) > 1Q 10 ^ 10 llQ eV _ (5 21) 

In the simple case when M*( mat ' and M* are of the same order this translates into the 
lower bound on M*. Note that this is several orders of magnitude below the upper limit 



( 5.17 ) imposed by the absence of strong coupling. 



It has been argued that considerably stronger constraints on the scale M* a come 



from the physics of ultra- high-energy cosmic rays (UHECR) Eq, 47, 48, E9]. These con 



straints are, however, less robust than ( 5.21| ) as they are sensitive to a number of assump- 



tions about the sources of cosmic rays and the UHECR chemical composition (e.g. the 
constraints essentially disappear if UHECR of the highest energies consist of heavy nuclei 
|49(1 ). Thus we do not discuss these constraints in the present paper. 

Finally, violation of Lorentz symmetry in the matter sector by operators of dimensions 
3 and four 4, such as the first and the second operators in ( 5.19| ), would lead to sizable effects 



even at low energies . On the other hand, the experimental constraints on these effects are 

29 Terms with odd powers of momentum can be forbidden by imposing the CPT invariance. In partic- 
ular, this allows to avoid bounds from the absence of vacuum birefringence discussed in M, B. It is also 
worth mentioning that modification of dispersion relation of the matter fields at high momenta is strongly 
suppressed in supersymmetric Lorentz-violating extensions of the Standard Model f42, [43]. 

It is worth stressing that without additional assumptions the khronon - matter couplings generically 
would be non-universal, i.e. species-dependent, implying, in particular, violation of the weak equivalence 
principle. 
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extremely tight || . One arrives to the conclusion that breaking of Lorentz symmetry in the 
matter sector at the level of dimension 3 and 4 operators must be highly suppressed. For 
some of the couplings this may be achieved by imposing discrete symmetries. For example, 



the first operator listed in (5.1£) can be forbidden by requiring CPT invariance. However, 
a stronger mechanism is needed to suppress all dimension 3 and 4 operators. We will 
mention one possibility in the Discussion section. For the moment let us just assume that 
such a mechanism exists. Then to the leading approximation the coupling of the khronon 
field to matter at the lowest-derivative level must be encoded in a universal effective metric 



9pLV 9[1V PU^Uy , 



(5.22) 



where j3 is a dimensionless constant. Clearly, this kind of coupling preserves Lorentz 
invariance of the matter sector. We concentrate on this type of coupling in what follows. 



5.3 Velocity-dependent forces and instantaneous interaction 

Before engaging into a systematic analysis of the observational constraints on the univer- 
sally coupled khronon let us make some preliminary estimates of the type of effects induced 
by the coupling ( |5.22| ). When /3 is small one expands the interaction terms and obtains to 
leading order in /3 



S v 



-mat 



I I ''" 



-g Ufj,u u 



rpflU 



where T^ u is the matter energy-momentum tensor. Expanding the khronon field as usual, 
<p = t + Xt one obtains the action for small perturbations including the source term 



/ 



A A x 



%a(9 lX f 



(\-l)(A X ) 2 )-f3xdiT 
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From this expression we read off the khronon exchange amplitude between two sources 
with energy-momentum tensors T^K and Tff!, 
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-lOi 



didj 
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a<9,j 
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T oj 

1)A ( 2 ) 



(5.23) 



This amplitude does not encapsulate all the new interactions due to the presence of the 
khronon field. In the systematic treatment Eq. ( |5.23j ) must be supplemented by the ampli- 
tudes involving khronon-graviton mixing and by the contributions coming from the mod- 
ification of the metric propagator due to the Lorentz symmetry breaking. This will be 
done in the next section within the formalism of the parameterized post-Newtonian (PPN) 



expansion. Here we note that the direct khronon exchange (5.23) dominates over other 
effects in the case 

a~|A-l|</3. (5.24) 



The simplicity of the amplitude (5.23) makes the analysis in this case much more trans- 



parent than in the full PPN treatment. Thus we first concentrate on the case (5.24). 

Let us consider explicitly the khronon-induced interaction for two point-like non- 



relativistic sources. 



In this case T?\ 
w 



m a v l a 5{^ — x a ), where x a , v a , a = 1,2 are the 
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coordinates and velocities of the two particles. In the non-relativistic regime one can neglect 
the time derivatives in the khronon propagator. After integration over spatial coordinates 
one obtains from ( |5.23 ) the following contribution into the two-particle Lagrangian, 



mim 2 /r vi • v 2 — (vi • r) (v 2 • ?) 



= ^_ ^ T± lj KZl VAIJ tl (5.25) 

x 8nMl A-l r ' v J 



where r = xi— x 2 and f = r/r. This interaction is quite peculiar: it depends simultaneously 
on the distance between the particles and on their velocities with respect to the preferred 
frame. Note that due the velocity dependence this term contributes non-trivially into the 
expressions for the conserved energy and momentum of the system. 



An interesting special case of Eq. (5.25) corresponds to the situation when one of the 

particles (say, particle 2) is much heavier than the other so that its velocity is approximately 

conserved (like, e.g., in the Sun - Earth system). One can show that in this case ( 5. 25] ) 

reduces to 

mim 2 (3 2 yj - (v 2 • f ) 2 

L * = 8^M| A=T r ' (5 ' 26) 

up to a total time-derivative. We observe that the dependence on the test particle velocity 
has dropped out. Equation (5.26) has a simple interpretation as a direction-dependent 



contribution into the gravitational potential of the heavy source. In the next section we 
will discuss the phenomenological constraints on this type of contributions within the PPN 
framework. 

An important remark is in order. The form of the khronon exchange amplitude ( 5.25| ) 



makes it manifest that the model involves instantaneous interaction. This is due to the fac- 
tor A -1 in the khronon propagator. This factor may cancel out in some special cases when 
the sources contain sufficient number of spatial derivatives, but not in general. Appearance 
of instantaneous interactions is a common feature of modified gravity models (more gener- 
ally, gauge theories) with broken Lorentz invariance [ll^, 50, [HI, 52 1. In a theory that does 



not aspire to a Lorentz invariant UV completion this does not pose an obvious obstruction. 
In particular, it certainly does not imply any problems with causality which is defined 
with respect to the preferred time slicing. Instantaneous interactions do introduce a kind 
of non-locality since far-away sources may affect the immediate future in any given local 
domain. However, the strength of the non-local effects seems to decay with distance and is 
suppressed both by post-Newtonian factors of v 2 and the (small) model parameters. Hence, 
it is not obvious if the presence of the instantaneous interaction leads to any significant 
constraints on the model. A more detailed investigation of related issues (the possibility 
to have any direct observational limits on instantaneous interactions, the implications for 
BH physics, etc.) is left for future. 

5.4 Universal coupling and Post-Newtonian Parameterization 

In the universally coupled case the effects of the khronon field are naturally interpreted as 
modification of the (universal) gravitational interaction between matter particles. The size 



of allowed effects is constrained by the existing tests of GR [37|. So we now turn to the 



bounds that the model has to satisfy in order to pass these tests. 
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The khronon sector is completely described by three parameters: (A — 1), a in the 
action ([|1]) and (3 in the effective metric ( |5.22| ) . It is convenient to make a field redefinition 
by trading the "Einstein-frame" metric g^v in favor of the effective metric g^ to which 
matter couples minimally. The result of this redefinition is readily obtained if we assume 
the parameters (A — 1), a, j3 to be much smaller than one; we will see that this assumption 
is justified by the phenomenological bounds which indeed require the above parameters to 
be small. Then to the leading order we write 

Using the identity 

one obtains the action 

S = -Mz- J d 4 xV=g{ i4) R + /SV^W + X'(V^) 2 + ouW/V^} , (5.27) 



where 



M' P 2 = Ml M + | ) , A' = A - 1 - (5.28) 



and we have omitted tildes over the new variables (we will always work with the redefined 
metric in the rest of this section so this will not lead to confusion) . Note that the appearance 
of the new parameter /3 in the khronon action ( 5.27| ) can be traced back to the presence 



of the parameter £ in front of the 3-dimensional scalar curvature in the unitary-gauge 
potential, see ( |2.6| ), (2.4). In the pure gravity theory this parameter does not have physical 



meaning as it can be eliminated, say, by a rescaling of the time coordinate, cf. the discussion 
after Eq. (|2.4|). This is no longer true in the presence of matter: as we are going to see 
physical observables depend on the value of f3. 

As already mentioned before, the action ( |5,27j ) has the same form as the action of 



the Einstein-aether theory [30| which has been extensively studied as a phenomenological 



model for violation of Lorentz invariance [31]. The difference of our model from Einstein- 



aether is that in our case the aether vector u^ is by construction hypersurface-orthogonal. 
As a consequence, in our case aether propagates a single longitudinal degree of freedom 
(khronon), while in general there are additional transverse modes. This implies that the 
results about Einstein-aether theory that are insensitive to the presence of transverse modes 
are also valid for our model 31 . 



31 The relation with Einstein-aether can also be used to derive the result of metric redefinition (5.22) 
in the general case when the parameters (A — 1), a, fi are not small. In Einstein-aether the substitution 
,22J) leads to the change of the coefficients in the aether Lagrangian; the corresponding formulas have 



been worked out in ]53j). The most general Einstein-aether action contains an additional term V^tivVw" 



compared to (3.27). In our case when the aether vector m m is hypersurface-orthogonal the four terms in the 



aether action are not independent [M implying that the term V^m^V^m" can be eliminated. In this way 
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The first set of constraints obtained in this way is related to the fact that the velocities 
of the modes propagated by the action ( |5.27| ) are in general different from 1, the maximal 
velocity of matter. Using the formulas for the Einstein-aether theory |3l| we obtain to first 
non-trivial order in parameters the velocity of helicity-2 modes (graviton) : 

c 2 g = 1 + $ , (5.29a) 

and that of the helicity-0 mode (khronon) : 

c 2 = >L±l . (5.29b) 



a 



Note that ( |5.29b| ) coincides with the velocity entering the dispersion relation ( |5.11 ). If the 



velocities of graviton or khronon are smaller than 1 relativistic matter particles will quickly 
loose their energy via vacuum Cherenkov radiation [54]. This is strongly contrained by the 



existence of high-energy cosmic rays. Thus we conclude that the graviton and khronon 
velocities must be larger or equal to 1 which yields the bounds 

/5>0, *±£>1. 

a 

Another constraint is obtained from the comparison of the gravitational constants 
appearing in the Newton law and the Friedmann equation governing the cosmological 
expansion. Again, the transverse aether modes do not play any role in these considerations, 
so we can directly apply the results from the Einstein-aether theory to our case. The 
Newton constant, which is defined as the coefficient in the Newton law for the gravitational 
force between two static masses, is related to the parameters appearing in the action ( 5.27| ) 
as follows Q: 

G N = gi . (5.30) 

8ttM' p 2 (1 - a/2) 

On the other hand, the cosmological expansion in the Einstein-aether theory is described 
by the standard Friedmann equation, 

H 2 = —G cosmP , (5.31) 

where H is the Hubble rate, p - the energy density of the Universe, but with a different 
proportionality coefficient 32 [31|, 



GcOSm "8vrM^(l + 3AV2 + /3/2)' (5 ' 32) 



one arrives at the action (5.27) with 



2 



M'/ = ^=, A' = (l-/?)(A-l)-/9 



which coincides with (5.25) for ^<1. 

32 The expressions (p.30|), (5.32) for the case — were derived directly from the action ([2.5]) in J25J 
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The discrepancy between Gn and G cosm is constrained by Big Bang nucleosynthesis [55] 



G 



1 



N 



< 0.13 



Barring accidental cancellations, this yields an order-of-magnitude bound on the parame- 
ters of the model, 

a, /3, X' < 0.1 . 

Stringent limits on any alternative theory of gravity come from the observational con- 
straints on the parameters of the parameterized post-Newtonian (PPN) formalism. Re- 
markably, in the Einstein-aether theory all PPN parameters except two are the same as in 
GR |5q| . The non-trivial parameters are called af PN , a PPN an d describe preferred frame 
effects related to breaking of Lorentz symmetry. We now argue that in the khrono-metric 
theory ( |5.27 ) all PPN parameters except af , a P are the same as in Einstein-aether. 



Indeed, besides a p , a P an d the parameters which vanish automatically in any theory 
described by a Lagrangian, there are three more PPN parameters: (3 , Y PN and £ PPN . 
The key point is that these three are determined from spherically symmetric solutions which 
are identical in the khrono-metric and Einstein-aether theories, see Appendix O Thus in 
the khrono-metric model at hand these parameters have the same values as in GR, 

p PPN = 1 PPN = l ^ ^PPN = Q _ 

One cannot use the results for Einstein- aether in the case of the parameters a PPN , 
a PPN : they describe effects related to the motion of the source with respect to the pre- 
ferred frame and are contaminated in Einstein-aether by the contributions of the transverse 
aether modes. These parameters can be defined as the coefficients in the linearized metric 



produced by a point source of mass m in its rest frame [57| 

m / (af PN - a PPN )v< 



hoo = —2Gn- 



r 



PPN 

a, „ m 



h oi = ^—Gn—v 1 
2 r 

VfX 

hij = —2Gn — Oij , 



af PJV (xV) 2 \ 
2 r 2 ) 


(5.33a) 




(5.33b) 




(5.33c) 



where r is the distance from the source and v l is the velocity of the source with respect 
to the preferred frame. Note that the contribution proportional to a PPN has the form of 



the direction-dependent gravitational potential encountered in Sec. 5.3. The current Solar 
system limits on these parameters are [J37f: 



\a PPN \ < lO" 4 , \a PPN \ < 10- 7 

I 1 I r^j I Z I r^- 1 



Derivation of a p , a P f° r the model ( 5.27| ) in the case a,/3,A' <C 1 is given in Ap- 
pendix ^; the result is 

„ppn A(„ offl „ppn (a - 2/3)(a - X' - 3(3) 

«l = - 4 (« - 2 P) 1 «2 = 2 r X , + o^ ■ I 5 - 34 ) 
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Note that both parameters vanish if a — 2/3 = 0. Another interesting case is (3 = 0, A' = a; 
this corresponds to the situation when the velocities of all modes in the theory are equal 
to 1, see Eqs. ( 5.29J ). In this case the parameter a.2 PN which is most tightly constrained 



vanishes. Barring these special cases and assuming a, (3, X' to be of the same order we 
obtain the bound 

a , , X' < 1(T 7 -r- 1(T 6 . (5.35) 

To the best of our knowledge, this is the strongest constraint on the parameters of the 
model that can be obtained at present. 

As discussed in Sec. 5.1 , the bound ( |5.35|) combined with the requirement that the 



theory is weakly coupled implies the upper limit fl5.17|) on the scale suppressing higher- 
order operators in the gravitational action. 

6. Summary and discussion 

In this paper we have investigated the self-consistency issues related to the scalar graviton 
modes in Hofava's approach to quantum gravity. We have considered three models of 
non-relativistic gravity differing by the symmetry group and the requirement or not of the 
projectability condition. In our study we extensively used the Stiickelberg formalism which 
makes the scalar modes explicit by encoding them into the khronon field: a scalar field 
with time-dependent VEV. This facilitates a lot the analysis of the scalar mode dynamics. 

In two of the considered models (the projectable version of the original proposal Q 
and a possible extension based on a smaller symmetry group), the scalar modes were 
found to exhibit pathological behavior, such as instability and strong coupling. These 
pathologies undermine the perturbative analysis. In particular, they invalidate the naive 
power-counting argument for renormalizability of these models. 

Several qualitative lessons can be extracted from these studies. While it is relatively 
easy to make the scalar modes well-behaved in the UV, this is much harder to achieve in 
the IR. The scalar gravitons tend to develop gradient instabilities which can be suppressed 
only by pushing the model parameters to extreme values. This, in turn, introduces strong 
coupling. Thus the primary goal in constructing a consistent non-relativistic gravity model 
is to stabilize the scalar modes. 

The analysis of Sec. ||] teaches us that reducing the symmetry to a smaller group than 
the FDiffs does not give any advantage in achieving this goal. Though this approach 
allows to improve the behavior of the scalar graviton of the projectable Hofava's model, 
it introduces yet other scalar modes which bring the pathologies back. Even relaxing 
the symmetry to the RFDiffs, that have as many local generators as the FDiffs, turns 
out to have the quite dramatic consequence of allowing for new operators that lead to 
tachyonic ghosts and therefore fast instabilities. Even if not present at tree level these 
operators would be generated by quantum effects. From this we conclude that the time- 
reparameterization symmetry contained in the FDiffs plays the quite important role of 
preventing the appearance of tachyonic ghosts in the model. 

Remarkably, the scalar mode is stable and also free from other pathologies in the 
third model which we analyzed. The latter is a natural extension of the original non- 
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projectable Hofava model obtained by including into the action all terms compatible with 
invariance under FDiffs and renormalizability by power counting. For appropriate choice of 
parameters the unique propagating scalar mode possesses stable dispersion relation in the 
entire range of spatial momenta. Moreover, the dispersion relation has nice properties both 
at low and high momenta. In the first case it is linear, u 2 oc p 2 , implying that the mode 
remains stable in any sufficiently smooth background 33 . On the other hand, the asymptotic 
form u! 2 oc p e of the dispersion at high momenta is compatible with the anisotropic scaling 
postulated in the UV. Therefore, the presence of this mode does not spoil the power- 
counting arguments which strongly suggest that the model is renormalizable. This implies, 
in particular, that strong coupling is avoided in the model by construction provided the 
coupling constants are chosen small enough. 

Encouraged by these results we studied some phenomenological aspects of the healthy 
model. Making use of the Stiickelberg formalism we have demonstrated that at low energies 
the theory reduces to GR interacting with an additional scalar field — the khronon. The 
time-dependent VEV of the khronon field breaks Lorentz invariance down to arbitrarily 
low energies. We observed that the structure of the low-energy limit of the model is similar 
to that of the Einstein-aether theory J30| , 31], even though the two theories are not exactly 



equivalent. The difference is due to transverse modes present in the Einstein-aether theory 
and absent in the healthy model of this paper. The transverse modes do not affect the 
form of homogeneous isotropic cosmological solutions and spherically symmetric solutions. 
This allowed us to directly apply to our case the bounds on Einstein-aether theory coming 
from the expansion history of the Universe, as well as to conclude that all but two PPN 
parameters in the healthy model coincide with those of GR. We have calculated the values 
of the remaining two PPN parameters, af , af , which characterize preferred frame 
effects. Current observational bounds on af , af significantly constrain the param- 
eter space of the healthy model (but do not rule it out). Combined with the requirement 
that the theory is weakly coupled these bounds translate into the upper bound on the scale 
of quantum gravity M* < 10 15 -r 10 16 GeV. Within the simplifying assumption that M* 
coincides with the scale of Lorentz symmetry breaking in the matter sector, astroparticle 
data constrain it from below: M* > 10 10 4- 10 11 GeV. [We stress though that this lower 
bound must be taken cautiously: its validity depends on the details of the matter sec- 
tor.] Thus the net result of our phenomenological study is that the healthy model can be 
tested experimentally using existing techniques, which we believe makes this model quite 
attractive. 

The above results suggest that the healthy model of this paper can serve as a starting 
point for construction of a viable renormalizable theory of gravity. Admittedly, a lot of 
open issues remain. Let us discuss some of them. 

The renormalizability of the model is yet to be demonstrated beyond the naive power 
counting. In principle, this amounts to an explicit analysis of the loop corrections to the 
action of the form fl2.5|) , (|2.6|). In practice though, this may turn out to be a formidable 
task, given the large (of order 100) number of allowed terms in the Lagrangian and the 

33 More precisely, in any background sufficiently close to Minkowski space-time equipped with the foliation 
by surfaces t — const. 
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complications related to the gauge symmetry of the action under FDiffs. A particularly 
subtle issue is a proper treatment of the time-reparameterization invariance. For consis- 
tency of the model this symmetry must be free of anomalies. Otherwise the symmetry 
would be reduced to that of the model of Sec. |], with the corresponding re-appearance 
of pathologies. An interesting development along these lines is the renormalization of the 
energy-momentum tensor of test fields in curved backgrounds in Hof ava-type theories J58| . 
It is shown that, in contrast to the relativistic case, this does not require counterterms with 
more than two time derivatives. 

Even if the model proves to be renormalizable, there will be a question if it is UV- 
complete or not. In other words, if it possesses a weakly coupled UV fixed point. An 
answer to this question requires the study of the renormalization group (RG) flow of the 
theory. Let us mention in this connection interesting recent results |59| about quantum 
electrodynamics (QED) in five space-time dimensions. It is demonstrated that 5d QED 
can be UV completed within Lorentz- violating framework by adding operators with higher 
spatial derivatives. The resulting theory is weakly coupled at all energies, possesses a 
weakly coupled UV fixed point with anisotropic scaling exponent z = 2 and flows to the 
usual 5d QED in IR (though, in general, with different velocities of photons and electrons). 

It will be interesting to assess the quantum properties of the healthy model beyond 
perturbation theory. A possible approach to this difficult problem would be to use the 
canonical formalism. It is worth noting that inclusion of the terms with derivatives of 



the lapse N into the potential (2.6) significantly improves the canonical structure of the 



theory compared to the original non-projectable version of the proposal. Indeed, due to 



these terms the lapse is no longer a Lagrange multiplier and the analysis [^, 24 1 unveiling 
the pathological structure of the constraints in the original Hofava model does not apply. 
Instead, in the model of this paper the Hamiltonian constraint obtained as the variation 
of the action with respect to the lapse, % = 4^ = 0, has non-vanishing Poisson bracket 
with the equation itn = 0, where 7Tjv is the canonically conjugate momentum for N. In 
other words, these two constraints form a second class pair and must be used to eliminate 
the variables N, ttn from the phase space [E3, pflj. Presence of second class constraints 
is an interesting feature of the model and its implications both for classical and quantum 
dynamics of the theory deserve a detailed study. 

The non-perturbative dynamics of the model can be also addressed semiclassically 
by developing perturbative expansion in the backgrounds of classical non-linear solutions. 
Clearly, the first step here is to find explicitly such solutions. Let us point out in this 
context the existence of cosmological solutions which are free from the initial singularities 34 



1 61, p2, 33 j. It will be interesting to understand if these solutions are stable. Another 
important issue is the structure of black hole solutions in the healthy model. Given that the 
high frequency modes in this model propagate with arbitrarily high velocities (moreover, as 



34 These solutions were obtained in the context of the original Hofava's proposal which suffer from the 
pathological behavior of the extra mode. However, the terms which were added to the original Lagrangian 
to obtain the healthy model vanish on spatially homogeneous configurations (they contain spatial gradients 
of the lapse). Thus the cosmological solutions of the original Hofava's model are also solutions of the healthy 
model of the present paper. 
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discussed in Sec. ^3, the theory involves instantaneous interactions), it is a priori unclear if 
the notion of black hole as a region of space surrounded by an event horizon makes sense. On 
the other hand, there certainly must exist some low-energy notion of horizon defined using 
low-frequency modes which have finite velocity. It will be interesting to understand the 
physical meaning of black hole entropy associated to this low-energy horizon. The situation 
is even more intriguing because generically, due to the violation of Lorentz invariance, 
different low-energy species will propagate with different velocities thus giving rise to a 



number of nested horizons. As discussed in [64, 65] one would expect this to lead to break 



down of black hole thermodynamics, which on general grounds implies violation of unitarity 
of the underlying theory. Thus it is vital for the proposal to understand if and how this 
paradox is resolved. It is worth mentioning that some spherically symmetric solutions of 



the healthy model have been reported in [66]. 



From the phenomenological perspective the major challenge for the non-relativistic 
gravity framework is a mechanism for emergence of Lorentz invariance in the matter sector 
at low energies. As discussed above this must happen with very high accuracy to satisfy 
existing experimental bounds. This poses a severe fine-tuning problem which is aggra- 
vated by the fact that Lorentz-violating parameters, such as differences of velocities of 



various matter species, run with the energy scale [p7| , 68|. A promising mechanism to avoid 
fine-tuning is related to supersymmetry. Indeed, it has been argued [42, ^3| that within 
Lorentz-violating extensions of the Minimal Supersymmetric Standard Model (MSSM) it 
is impossible to write any Lorentz-violating operator of dimension less than 5. Dimension 
5 Lorentz-breaking operators can be further forbidden by imposing discrete symmetries, 
e.g. CPT; in this case Lorentz breaking starts at dimension 6. In other words, given su- 
persymmetry, Lorentz invariance emerges as an accidental low-energy symmetry. Eventual 
supersymmetry breaking gives rise to Lorentz-violating effects, but these are suppressed 
by the small ratio of the soft supersymmetry breaking masses m so ft to the scale M* of 
the higher-order Lorentz-violating operators. For example, the coefficients of dimension-4 
Lorentz-violating operators generated from the operators of dimension 6 upon supersym- 
metry breaking are of order (m^t/M*) 2 . These are comfortably within the experimental 
bounds for m so f t ~ 1 TeV, M* ~ 10 15 GeV. Interestingly, supersymmetry also suppresses 
the contributions with higher powers of momentum in the dispersion relations of matter 



particles |42|, 43 1 which weakens significantly the lower limits on M* coming from astrophys- 
ical observations. Needless to say, realization of this scenario for the emergence of Lorentz 
invariance requires supersymmetrizing the non-relativistic gravitational action ( |2.5[) , (]2. 6]). 
The first step in this direction would be to super symmetrize the low-energy limit of this 
action, Eq. ( 5.27| ). This essentially amounts to constructing a supersymmetric action for 



the khronon sector as the gravitational part coincides with GR and can be made super- 
symmetric in the standard way. 

The present work makes only first steps in the study of observational consequences of 
the healthy model. Many topics which have not been touched in this paper deserve a thor- 
ough investigation. These include implications of the model for emission and propagation 
of gravitational waves, dynamics of binary pulsars, spectrum of CMB perturbations and 
structure formation. Interesting results about the dynamics of the cosmological perturba- 
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tions in the healthy model were reported recently in |39|, [7(]] . 
We plan to return to some of the above questions in future. 

Note added. When this paper was in preparation the article [J7l]] appeared which also 
uses the Stiickelberg formalism to study the properties of the healthy model. Where there 
is an overlap, our results agree. However, Ref. [71] pushes the analysis using the decoupling 



limit to the regime of high energies (higher than M*) and argues that the terms with the 
higher spatial derivatives do not resolve strong coupling of the scalar mode. This seems 
to contradict the conclusions of the present paper. In fact, pushing the decoupling limit 



to high energies, as done in JnJ, corresponds to going beyond its range of applicability. 
Indeed, the limit considered in |71[ corresponds to throwing away all the higher-derivative 
mixings between the khronon field and the metric. On the other hand, as admitted in the 
published version of the |7l|, these terms are crucial to ensure the correct UV behavior 
uj 1 oc p 6 of the dispersion relation of the scalar mode, required by the anisotropic scaling 



with z = 3 and satisfied in the full theory, see Eq. (5.5). Instead, for the decoupling limit 



considered in [[n]] the UV asymptotic of the dispersion relation for the khronon field is 
to 2 = const. Thus the loop integrals over spatial momenta instead of being suppressed 
become even more divergent than in the relativistic case. Clearly, this only represents a 
failure of the limit considered in |7l]] and does not imply the inconsistency of the theory. 
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A. Proof of the no ghost theorem 



In this appendix we prove the proposition formulated at the end of Sec. 2.2. For simplicity, 
let us first consider the case when the Lagrangian depends only on the normal vector u M 
and its first derivatives. These derivatives enter into the Lagrangian through the extrinsic 
curvature Ku, v and the acceleration vector a M , 

This case covers all terms in the general Hofava-type Lagrangian ( [2.5D except those involv- 
ing spatial derivatives of the 3-dimensional tensor Rij or of the acceleration a,. 
We start by expanding u M up to quadratic order, 
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where 

^\ = ^=P;d vX , (A.l) 

6&\ = -lu,{P» x d xX f - jP'dxX u"d vX ■ (A.2) 

Here bar refers to the background values of the fields. The crucial observation is that in 
the background (|2.20| ) the first variation ( A.l ) of the normal vector does not contain time 



derivatives of Xi while the second variation ( |A.2 ) contains only first time derivative. The 



next step is to consider variations of the extrinsic curvature and acceleration 
S^K^ = S^P X V x u u + P x V x S^u u , 

8^K, V = 26Mpi V x 5^u u + 6® pA Vx u v + P x V A ^V , 
5^a p = 25^u x V x 5^u p + 5^u x V A u M + u x V x 5^u^ . 



One observes that in the background ( 2.20J ) 5^ l 'K,^ u does not contain time derivatives of 



X', S^'ap, 5^ 2 'K,^ V contain one time derivative, and d^'a^ - two time derivatives. The 
quadratic Lagrangian consists of two types of terms: the terms containing the first vari- 
ations of the fields squared, and the terms linear in the second variations. For example, 
writing down explicitly the terms with the variation of the acceleration we obtain 



(2) _ 1 d 2 C (1) (1) ld£ (2) 

~ 2 da da M 2 da M ' ' ' ' 

Clearly, both these terms contain at most two time derivatives of x- I n the case when (p 
does not satisfy equations of motion one has to consider also the linear variation of the 
Lagrangian D- 1 '. Using the same reasoning as before one concludes that C^ l > contains at 
most first time derivative of x- This completes the proof. 

It is straightforward to generalize the above proof to include the dependence of the La- 
grangian on higher spatial derivatives. In covariant language a spatial derivative translates 
into the operator P o V o P, which is purely spatial (i.e. does not introduce further time 
derivatives). For example, the object VjOj takes the covariant form P x V x (Pva p ). Using 
the same reasoning as before one can show that the first (second) variation of this type of 
objects contain at most one (two) time derivatives of %. 

Note that the above proof uses in an essential way the invariance of the action under 



the transformations (2.15): this forces the Lagrangian to depend on the khronon field only 



via Un and its derivatives. Recall that this invariance stems from the invariance of the 



original action fl2,5|) under FDiffs. As a consequence the proof also goes through for the 



projectable model (|2.2j ) which obeys this symmetry. Indeed, as pointed out in Sec. ^1 
the projectable model is recovered from the theory (|2.5|) , ( |2.6|) by taking the limit a — > 00. 
Thus the order of time derivatives in the covariant equations for these two theories coincide. 
On the other hand, the proof does not apply to models with reduced symmetry, such 
as the theory fl2.7|) , because in these cases the Lagrangian contains additional dependence 
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on the tp gradients. We will see in Sec. |4| that this makes the equation of motion for the 
khronon fourth-order in time; as one can anticipate, this will lead to certain pathologies of 
the theory. 

B. Two faces of the tachyonic ghost 



The purpose of this appendix is to clarify the following puzzle. In Sec. B we have shown 
using the Stiickelberg formalism that the RFDiff model includes in its spectrum a tachyonic 
ghost. On the other hand, in the unitary gauge the action (|2.7| ) of the model contains only 
two time derivatives and according to the standard lore one would not expect any ghosts 
in this picture. This seems to contradict the results of the Stiickelberg analysis. We are 
going to show that the contradiction is removed when one properly formulates the physical 
questions to assess the effects of the tachyonic ghost. 

We start by writing the unitary gauge Lagrangian for the sector of scalar perturbations 
around Minkowski background. Substituting the decomposition ( ^S| ) in the action ( p.7| ) 
and integrating out the non-dynamical fields B and E one obtains 



£ 
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III 
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2-0AV> + 44>Ai/j 
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2 M l 
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* (Atf + % A^ + A(A^ + &** + |>3, + * M «, 



(B.l) 



where M\, M A1 , M a are defined in ( 2.1C| ) and M 4 is the coefficient appearing in the 
expansion of the potential V(N) in ( |2.8[ ) to quadratic order in eft; these parameters have the 
same meaning as in the Stiickelberg action (fDj) . In deriving ( |B.1| ) we assumed for simplicity 
M\ <C Mp (this is the case relevant for comparison with the Stiickelberg analysis). Finally, 
the constants f n , g n are related to the coefficients of the higher-derivative terms in the 
potential (|z 



The Lagrangian (BA) clearly describes two propagating degrees of freedom 
which matches with the Stiickelberg analysis of Sec. ||. What seems to be different is that 
the kinetic energy of both modes can be made positive by choosing M\, M A1 > 0, while in 
the khronon language one of the modes is a ghost. This said, let us proceed and find the 
dispersion relation for the modes. Neglecting the higher-derivative terms we obtain, 
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(B.2) 



Clearly one of the modes exhibits gradient instability. In the decoupling limit Mp — > oo 
the dispersion relations ( B.2 ) coincide with the expressions (|4.6|) of Sec. |I|. Thus at the 
level of the dispersion relations the unitary gauge and the Stiickelberg descriptions match. 
Let us now ask if there is any physical setup where one could distinguish between a 
mode with gradient instability (simple tachyon) and a mode which besides the gradient 
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instability also has a negative kinetic term (tachyonic ghost). The standard definition of 
the ghost as "a field whose Hamiltonian is not positive definite" is not very useful: the 
Hamiltonian is not sign-definite in both cases. Moreover, the two cases can be related by a 
canonical transformation. To illustrate this point consider a toy Lagrangian representing 
a single Fourier mode of a tachyonic ghost, 

^2 ^,2^,2 

(B.3) 



L'ghost — 


if p 2 7] 2 
2 2 


The corresponding Hamiltonian reads 




■" ghost 


ir 2 p 2 rj 2 
2 2 



where n is the canonically conjugate momentum for rj. The canonical transformation 

7r = prj , f\ = ir/p 
casts this into the Hamiltonian of a simple tachyon, 

_ 7f 2 pV 

■" tachyon „ ~ • 

One may still try to distinguish if a mode is ghost or not by the sign of the residue at 
the pole in the one-particle exchange amplitude considered as function of w 2 ; this must be 
negative for the mode to qualify as a ghost. However, in the case of the tachyonic ghost 
the sign depends on the type of the coupling to the source used to define the amplitude. 
Taking again the toy model (p.3|) as an example consider two couplings: 



^soiree = V^i and L^ 2 J urce = 77E2 • 



•Asi oc Ei — — — - Ei 



In the first case the one-particle exchange amplitude reads 

-1 

~P l 
and the residue is negative. However, in the second case 

Ay, 2 oc £ 2 —5— — n E2 
u z + p A 

and the residue becomes positive: — uj 2 = p 2 > 0. Note that this ambiguity is related to 
the instability of the mode; it is absent for the case of a ghost with a stable dispersion 
relation when the pole lies at positive co 2 . Thus we conclude that for a mode with gradient 
instability there is no unambiguous way to tell if it is a ghost or not. Rather, this notion 
makes sense only for a given coupling of the mode to the source. 

It is instructive to trace explicitly the agreement of one-khronon exchange amplitudes 
calculated in the Stiickelberg and unitary gauge pictures. To this aim we need to specify 
the source. This is easier to do on the Stiickelberg side where we couple the khronon to a 
scalar field E. The khronon field must enter with derivatives so we write, 



S„ 



I tfxJ^ggrdpvdvZ. (B.4) 



40 



In terms of the khronon perturbations this takes the form 

Ssource = / d A X (~X + A*) E 



From this expression and the quadratic action (4.5) one reads off the khronon exchange 
amplitude 

I 2 _ 2n2 

Aj, oc S * ^- E , (B.5) 



where 



L> = M^w 4 + M 2 w 2 p 2 - M|p 4 + 2MV . 



One observes that this amplitude is a sum of two contributions, with positive and negative 
residues at the poles. In this sense one of the khronon modes is indeed a ghost. 

Let us now see how the same amplitude is recovered in the unitary gauge. Fixing ip = t 



and substituting the decomposition (23) into the source term flB.4|) we obtain 



f d 4 x (-<j) + 2iP + E + VAB) E . 



The next step is to integrate out the non-dynamical fields B, E. Importantly, this produces 
a contribution into the Lagrangian which is quadratic in E. Omitting the higher-derivative 
terms the resulting Lagrangian reads 

A/f 2 ■ A/f 2 A/f 2 ■ V 2 

4?/ = ^ + *4M + ^{d l(t >? + -ftf + MV - {4> + tf)E + Aj , 

A 

where for simplicity we have taken the limit Mp — > oo (this corresponds to the decoupling 
limit in the Sttickelberg picture) and introduced the rescaled variable 

4, = (2M 2 P /M 2 X ) V . 

This Lagrangian leads to the following propagators: 

Finally, the khronon exchange amplitude reads, 

„4 E ocE (w 2 ((<^)+2(<^) + (^))-M^ 2 ) E. 



Combining everything together one obtains the result ( B.5 ). Thus we find that the unitary 
gauge calculation leads (in a quite non-trivial way) to the same result as in the Stiickelberg 
formalism, as, of course, it should be. Technically, in the unitary gauge the ghost pole 
appears due to the structure of the coupling to the source. 
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C. Stability bounds for the healthy model 

Here we present conditions on the parameters of the healthy model of Sec. ^ imposed by 
requiring the stability of the scalar mode. We formulate them in terms of the coefficients 
a ) fn, 9n appearing in the quadratic Lagrangian ( |5.1| ). 

The stability requirement is expressed by Eq. (|5.6|) with the polynomials P{x) and 



Q(x) defined in ( |5.3| ), (5.4). This implies the following necessary conditions: 

g\ - gm >o,2-a>o, 

g 3 > , a > , h < 2V®93- 

Deriving the full set of necessary and sufficient conditions is quite cumbersome. Instead 
we provide two different sets of sufficient conditions. The first possibility is requiring that 



all the monomials in P(x) are positive definite. Apart from the previous conditions ( C.l ) 
this yields the constraints 

51/3 + 53/1 -252/2 > , 

/ 2 2 - %2 - /1/3 - 2g 3 -a 9l >0, (C.2) 

2/3 + «/i + 4/ 2 > , 
Another option is to write P(x) as 

P(x) = {{gl - gm) 1/4 x + (4 - 2a) 1 / 4 ) 4 - x(c 2 x 2 - Cl x + c ] 



and require the quadratic polynomial inside the last bracket to be positive at x < 0. This 
translates into the constraints 



where 



c 2 > , ci > -2^/qJc2~ , c > , (C.3) 

c 2 = fm + fm - 2f 2 g 2 + 4( 52 2 " 5i5s) 3/4 (4 - 2a) 1 / 4 , 

ci = /| - 4g 2 - A/3 - 2g 3 - 9l a - 6(g 2 2 - gm) 1/2 (4 - 2a) 1 / 2 , 

c = 2/ 3 + 4/ 2 + ha + \{g\ - 5 i5 3 ) 1/4 (4 - 2a) 3 / 4 . 



Note that the two sets of bounds ( |C.2|) and (|C.3|) overlap but none of them contains the 
other. 

Clearly, the necessary conditions ( |C.1| ) can be complemented with any of the sufficient 
conditions ( |C.2| ) and ( |C.3j ). To demonstrate that the parameter space restricted by the 
stability bounds is not empty let us give an explicit example. It is straightforward to verify 
that the set of parameters 

a = #3 = 1-5 , /i = f 2 = ~h = ~g\ = -52 = 2 



satisfies the constraints ( |C.l| ), ( C.2 ) and thus leads to stable dispersion relation of the 
scalar mode. 

The bounds presented above can be translated directly in terms of the parameters Ai, 



E>i, Ci and Di in the original potential (2.6) but we do not do it here 
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D. Spherically symmetric solutions in Einstein-aether and khrono-metric 
theories 

In this appendix we demonstrate that spherically symmetric solutions of the khrono-metric 
model ( |5.27D are identical to those of the Einstein-aether theory. Let us consider the 
equation of motion for the khronon field coming from varying the action ( |5.27 ) with respect 
to the field x 

, (D.l) 



V^J" 



where 



J' 1 



1 5S 



At the same time the equation of motion for the aether is obtained by varying ( |5.27 ) with 
respect to the field u v and reads 



V 



1 SS 



0. 



(D.2) 



In deriving this equation one has to take into account the constraint u^u^ = 1: it leads 
to the appearance of the projector P„ on the l.h.s. Finally, the energy-momentum tensor 
appearing in the Einstein equations for the khrono-metric theory can be written as 



T, 



SS 



fJ,U 



-9 \ $9 



fJ,V 



1 SS 

2S~^' 



where the second term comes from the explicit dependence of the vector Ua in the khronon 



theory on the metric, see Eqs.(2.13), (2.14). This coincides with the energy-momentum 
tensor of the aether. To obtain the second term in this case one again has to take into 
account the constraint u^u^ = 1 p4f . 

Any spherically symmetric configuration of aether is automatically hyper surf ace-or- 
thogonal implying that any spherically symmetric solution of the Einstein-aether theory is 
a solution for the khrono-metric theories [34|. The converse is less obvious as the equation 
of motion ([D.lD of the khronon field contains an additional derivative compared to the 
aether equation ( p.2|) and thus, a priori, admits more solutions. However, for spherically 
symmetric configurations ( D.l ) implies ( p.2[) . Indeed, the current J^ obeys the relation 



Ufl J» = . 



Hence in the unitary gauge (2.12) its component J° identically vanishes imlying that the 
corresponding charge 

Q= fd 3 x^J° 

is identically zero. On the other hand, the time derivative of Q is equal to the flux of 
the spatial component J 1 of the current through the 2-sphere at spatial infinity 35 . In 



35 We assume that the 3d surfaces ip = const do not have holes. This is the case if these surfaces form a 
regular foliation of the whole space-time. 
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spherically symmetric situation this implies that the current itself is zero which brings us 
to the equation ( p,2| ). Combining this result with the equality of the energy- momentum 
tensors we conclude that spherically symmetric solutions in the khronon and Einstein- 
aether theories are indeed identical. 



E. PPN parameters af , a? f° r the healthy model 

In this appendix we derive the formulas for the PPN parameters a PPN , a PPN in the 
khronon theory (|5.27). We assume the couplings a,/3, A' to be small and perform calcula- 



tions to the leading order in these couplings. We consider the metric produced by a point 
source of mass m in its rest frame. This frame does not coincide with the frame defined by 
the preferred foliation. Hence the background value <p of the khronon field in this frame 
differs from the coordinate time. Using the reparameterization symmetry (|2.15| ) we fix 



r 



Vl + v 2 t + v i x i , 



where v % is the velocity of the source with respect to the preferred frame. This corresponds 
to the background value of the vector u M 



u = V 1 + v 2 , Ui = v l . 
The source perturbs the metric and the khronon. One writes, 

9fiu = rjtiv + hyv , Lp = (p + x ■ 
To the leading order hu,v is given by the standard Newtonian expressions, 

h<$=26(r), 4?=0> h$=2<Kr)6 ij , (E.l) 

where 

Sir) = =j- . (E.2) 



Our goal is to find corrections to ( [E.lD in powers of v. 
It is convenient to introduce the following notations 

$ = haPifhiv, , v p = vfP^h^ , h Xp = P^PpKv , (e.3) 

a" = u^ , d^ = P»d u , D x = d^d ±u . 

The indices here are raised and lowered using the Minkowski metric r] pv . Expanding to 
linear order in perturbations we obtain 

u^ = Up + d p x + u^ , 
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Substituting these expressions into the action for the khronon sector (last three terms in 
( 5.27Q ) we obtain at the quadratic level, 

$x = "^ |d 4 x|/3((D^) 2 - 2D ± X a^V - dhd^H^ + \{d^f 






(E.4) 



where H = H" . This yields the equation for the khronon perturbation x, 

(E.5) 
Variation of the action ( |E.4| ) with respect to the metric perturbation hp U gives linearized 
khronon energy- momentum tensor: 

~6S V 



rppLV 






(uVPZ + WP^ 6 ^ 



A 



2 p»pv 5S X- 



(E.6) 



5$ ' ' A/ 5V X ' x p 8H Xp 

where we have used the decomposition 

hyu, = 2u^u u ® + u,j,V v + UvVfj, + Hfj, v . 

Evaluating the variations entering into ( |E.6| ) we obtain 
SS V 



5$ 



Mj?a(0 ll n- L x-n- L $ 



5S- 



x 



SV X 



-M 



/2 



(A' + p)d^ x ~ (A' + d ±X d^V 



(3 



A' i x. 



P, 



5S y 



5H 



Xp 



-Mj 2 



2 2 2 p 

y 

2 



J#(q\U l x ~ d^V* + -{d^) 2 H) 



P. 



P, 



_ P d \\ d ^ d ^ + P d \\ d ± x V P + !i d \\ d ^v x - ^(d\\?H x <> 

2 4 4 4 



(E.7a) 



(E.7b) 



(E.7c) 



The rest of the calculation proceeds as follows. One inserts the Newtonian metrics 



K|) into Eqs. (j|J) and find the potentials $ (0) , V^ j and H ( ^J . The latter act as 



x P 



source for the khronon perturbation x i n Eq. (|E.5| ). At the next step one combines the 
khronon perturbation found from ( |E.5[) and the Newtonian expressions for the potentials 
into the khronon energy-momentum tensor ( |E,6| ). This tensor substituted in the r.h.s. of 
the Einstein's equations determines the correction to the metric: 



Ah™ - — , , 

^ v ~ M' 2 \ x '"' 



2 r lt iV x 



(Ei 
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Here we have imposed the harmonic gauge, 

dM lu - -d u h = 

and have used the fact that the metric is static. Note that the first order in the post- 
Newtonian approximation requires to find h 00 , h 0i and h %3 components of the metric with 
the accuracy 0(v 2 ), 0(v) and 0(1) respectively. This implies that we need to determine 
T£ v and T x to order 0(v 2 ), T® 1 - to order 0(v), and T% - to order 0(1). 
Expanding up to terms 0(v 2 ) we obtain: 





uq = u° = 1 + v 2 /2 , m = -u l = v l , 


p° 


= -v 2 , p* = -po = v i ) pj = s l j + v l v j 


9" 


= -v l 3i , d£ = v i 8 l , di = di + vVdj 




D 1 - = -A - v l v j d t dj , 



where we have used that derivatives act on static configurations. Substituting these ex- 
pressions together with ( jE.l|) into flE.3|) we find 



$(0) = (1 + 2v 2 )cj)(r) , 

F o (0) = -4i> 2 ^(r) , T// 0) = -4«V(r) , 

#$ = 2v 2 <j>(r) , fl£> = 2uV(r) , H\f = (2<% + 6vV)0(r) , 

F (0) = (-6 - 4u 2 )0(r) 

The khronon equation flE.5| ) takes the form, 

(A' + /3)A 2 X = (a - V - SpydiA^r) , (E.9) 

where on the l.h.s we have neglected terms 0(v 2 A 2 x) as they are of higher order in v. 
From ( |E.9| ) we find 

a- A' -3/3 j . . 
Ax = A/ ~ -g ^9^(r) , 

meaning that x is °f order 0(i>). 



Let us estimate the orders of the variations (E.7). One finds by inspection that -?# is 
of order O(l), -^jr- - at most of order 0(v), jjp — at most of order 0(v 2 ). This allows to 
simplify the khronon energy-momentum tensor. To the required orders we have: 

T xW = -(l + v^ + 2v' S -^, (E.lOa) 

T *—'& +"*$;■ (E - 10b) 

r *°-= -"'§+§• < E - 10c > 

T xij = 0. (E.lOd) 
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One evaluates the variations appearing in these formulas, 



s_s 2L _ 

5S X __ 
d%3 8H,- 



-M$a 



( I + 2v 2 )A(f ) (r) + ^-^-vVdidj^r) 



-M'£ 



-M 



2Pv i A(f)(r) + (a - 20)it didrfir) 

, 2 (a-2/3)(3A' + /3) 



'j 



2(A' + /3) 



-v l v 3 didj(j)(r) . 



(E.lla) 
(E.llb) 
(E.llc) 



Inserting ( |E.11| ) into ( |E.10 ) and substituting the result into ( |E.8[) we find the equations for 
the first order corrections to the metric, 

A/4o = (« + 4(a - 2/3)t> 2 ) A0(r) + ( a ~ 2 ^~ q X ~ 3/ \ v<9^-(/>(r) , (E.12a) 



A' + /3 



,W 



A/i^ ; = 2(a - 2 / 9)t; 1 A0(r) - 2(a - 2/3)v ] did j( f>(r) 
Ah[f = aSijA^r) . 



(E.12b) 
(E.12c) 



It is straightforward to solve these equations for the explicit form flE.2|) of the function 
4>(r). Note that the second term on the r.h.s. of ( E.12b| ) can be removed by a time- 
independent gauge transformation. Combining the result with the Newtonian expressions 
( |E.1| ) we obtain the metric ( 5.33J ) with 



Gn 



1 



8irMg 



' + ! 



(E.13) 



and the PPN parameters af , af quoted in ( 5.34| ). The expression ( |E.13| ) coincides 
with (5.30) to linear order in a. 
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